The  Four-Dimensional  Fire-Induced 
Transmittance  and  Turbulence 
Effects  Model:  FITTE 

by  Dorothy  Bruce 

Battlefield  Environment  Directorate 


ARL-TR-273-6 


August  1995 


19960327  051 

Approved  for  public  release;  distribution  is  unlimited. 


me  QUALTIY  IMBPECl’SD  X 


NOTICES 


Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of 
the  Amiy  position,  unless  so  designated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in  this  report  is  not  to 
be  construed  as  official  Government  indorsement  or  approval  of  commercial 
products  or  services  referenced  herein. 


Destruction  Notice 

When  this  document  is  no  longer  needed,  destroy  it  by  any  method  that  will 
prevent  disclosure  of  its  contents  or  reconstruction  of  the  document. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704^188 

Public  reporting  burden  for  this  collection  of  information  Is  estimated  to  average  i  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gattv^ng  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
colleaion  of  information,  including  suggestions  for  r^ucing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  121 S  Jefferson 
Davis  Highway,  Suite  1204.  Arlington,  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

August  1995  Final 

4.  TITLE  AND  SUBTITLE 

The  Four-Dimensional  Fire-Induced  Transmittance  and  Turbulence 
Effects  Module:  FTITE 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Dorothy  Bruce 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

U.S.  Army  Research  Laboratory 

Battlefield  Environment  Directorate 

Attn:  AMSRL-BE-S 

White  Sands  Missile  Ranee.  NM  88002-5501 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ARL-TR-273-6 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 

U.S.  Army  Research  Laboratory 

2800  Powder  Mill  Road 

Adelphi,MD  20783-1145 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

ARL-TR-273-6 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

A 

13.  ABSTRACT  (Maximum  200  words) 

The  four-dimensional  fire-induced  transmittance  and  turbulence  effects  module  (FITTE)  calculates 
transmittance  and  laser  beam  spread  and  wander  for  propagation  paths  through  battlefield  vehicular 
or  hydrocarbon  fires  and  also  calculates  radiance  from  the  fire  and/or  fire  plume.  It  can  perform  the 
transmittance  and  radiance  calculations  for  a  rectangular  array  of  points  and  create  an  output  file 
which  can  be  used  to  add  images  of  battlefield  fires  into  scenes. 

This  document  contains  both  a  technical  description  of  the  model  and  a  user’s  guide. 

14.  SUBJECT  TERMS 

fire,  transmittance,  turbulence  effects,  radiance,  beam  wander,  beam  spread, 
emissive  sources 

SBS 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

— Unclassified  _ 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  Z39-18 
298-102 


Acknowledgments 

The  author  wishes  to  acknowledge  helpful  discussions  at  the  start  of  this  work 
with  Dr.  Donald  Walters  and  Dr.  Kenneth  Kunkel  about  atmospheric  turbulence 
theory  and  measurements,  and  recent  discussions  with  Dr.  Donald  Hoock  about 
fractal  representations  of  fluctuational  quantities.  She  is  also  indebted  to  Dr. 
Charles  Bruce  for  data  and  discussions  of  their  interpretation,  and  to  Wendell 
Watkins  and  his  coworkers  for  transmittance  and  calibrated  imagery  data. 
Members  of  the  NATO  Research  Study  Group  15  and  their  technical  support 
associates  provided  many  interesting  discussions  and  extensive  data  that  were 
useful  in  model  evaluation.  Special  thanks  are  extended  to  Dr.  Dieter  Clement 
of  the  Forschungsinstitut  fiir  Optik  for  arranging  to  share  the  model 
enhancements  they  supported. 


1 


Contents 


Acknowledgments .  1 

1.  Introduction  .  7 

1. 1  Names  of  Things .  7 

1.2  Availability  .  8 

1.2.1  Mailing  Address .  8 

1.2.2  Phone  and  Electronic  Mail  .  . .  8 

2.  Background .  9 

2. 1  Conceptual  Overview  of  FITTE  .  11 

2.2  The  Physical  Models .  11 

2.2.1  The  Fire  and  Fire  Plume  Models .  11 

2.2.2  The  Radiative  Transfer  Model .  25 

2.2.3  Fire  Plume  Turbulence  Effects  on  Laser  Beam  Propagation .  29 

2.2  Model  Implementation . 34 

2.3.1  FITTE .  34 

2.3.2  FGLOW:  Fire/Plume  Radiant  Emission  and  Transmittance .  38 

3.  Caveats .  41 

3.1  Grade  of  Software  .  41 

3.2  Model  Failure  .  41 

3.3  Verification  Tests  .  42 

3.3.1  Fire  and  Plume  Properties .  42 

3.3.2  Aerosol  Optical  Properties .  44 

3.3.3  Plume  Distribution  Functions .  46 

3.3.4  Effects  of  Thermal  Plume  Turbulence  on  Laser  Beam  Propagation  .  49 

4.  Operations  Guide .  51 

4.1  Input .  51 

4.2  Output .  61 

4.2.1  Parameters  Returned  to  the  Calling  Routine .  61 

4.2.2  Output  Printed  or  Stored  in  Files .  61 

5.  Sample  Runs .  63 

5.1  Example  1:  A  Scenario  1  Calculation  for  a  Single  Wavelength .  64 

5.2  Example  2:  A  Scenario  2  Calculation  for  a  Single  Wavelength .  65 

5.3  Example  3:  A  Scenario  1  Calculation  for  a  (Visible)  Waveband .  66 


3 


5.4  Example  4:  A  Scenario  1  Time-Series  Calculation  for  a  Single 

Wavelength  .  67 

5. 5  Example  5:  A  Scenario  1  Calculation  With  Variation  of  Fire 
Parameters  and  Printout  of  Spectrally  Resolved  Results  for  a 

Selected  Instrument  Response  Function .  69 

5.6  Example  6:  A  Single  Wavelength  Scenario  3  Calculation .  71 

5. 7  Example  7:  A  Scenario  3  Waveband  Calculation .  73 

5.8  Example  8:  A  Time-Series  of  Single  Wavelength  Scenario  3 

Calculations . 75 

5.9  Example  9:  Internal  Cycling  With  Variation  of  Wind  Direction  to 

Produce  Five  Output  Files .  77 

References  .  g5 

Acronyms  and  Abbreviations  .  89 

Appendices 

Appendix  A.  Fluctuation  Time-History  Data  Files .  91 

Appendix  B.  Auxiliary  Programs  for  Visualization  of  FGLOW  Output .  95 

Appendix  C.  Modification  of  Fire  Parameters  .  99 

Appendix  D.  Model  Organization .  105 

Distribution  .  II3 


Figures 

1 .  Side  view  contour  plot  of  a  fire  plume  in  the  model  coordinate  system. 


The  plume  edges  shown  are  at  the  e'^  points  of  a  conical 

Gaussian  distribution.  The  centerline  is  sketched .  15 

2.  Plume  distribution  function  behavior  as  a  function  of  distance  along 

the  centerline  from  the  fire  .  18 

3.  Equivalent  top  hat  and  Gaussian  distributions  .  19 

4.  Contour  plot  of  a  space  and  time  varying  (4D)  plume.  Tic  marks 

show  5-m  spacing.  The  dashed  line  is  the  sketched  centerline .  26 

5.  Geometry  for  the  single  LOS  scenarios.  The  helicopter  carries 
a  laser  designator  for  scenario  1  or  a  laser-guided  munition  for 

scenario  2 .  35 

6.  Graphs  of  model  output  data  for  a  time-series  calculation  at  10.6  pm. 

The  dashed  lines  show  the  values  predicted  by  the  mean- value  model  .  36 


4 


7.  Top  view  of  scenario  3  (FGLOW)  geometiy.  Situation  shown  is  for 

the  240°  wind  direction  of  example  9 .  39 

8.  Comparison  of  measured  and  modeled  plume  elevation  angles  .  44 

9.  Comparison  of  measured  and  modeled  extinction  coefficients .  45 

10a.  Transmittance  plots  created  with  data  from  the  example  8  simulation. 

Wavelength  is  10.6  pm.  Tic  spacing  is  5  . .  76 

10b.  Radiance  plots  created  with  data  from  the  example  8  simulation. 

Wavelength  is  10.6  pm.  Tic  spacing  is  5  . .  76 

1 1 .  Apparent  temperature  plots  for  various  wind  directions. 

Prepared  from  the  data  files  of  the  example  9  simulation. 

Wavelength  is  3.6  pm .  ^3 

Appendix  Figures 

B-1.  Greyscale  print  of  plume  optical  density  produced  with  XS2PRT  .  98 

D-1.  Diagram  of  FITTE  model  structure.  Dashed  lines  indicate  conditional 
calls.  Functions  called,  except  EXPOF,  are  listed  in  parentheses  in 

subroutine  boxes . 1^7 


Tables 

1.  Values  defining  model  fires .  12 

2.  Particulate  optical  parameters .  23 

3.  Comparison  of  albedo  values .  46 

4.  The  REFD  card .  32 

5.  The  SRCL  card .  53 

6.  The  SCEN  card .  53 

7.  The  METD  card  .  54 

8.  The  DETD  card  . 54 

9.  The  TARG  card  .  54 

10.  The  MOLS  card  .  55 

11.  The  BAND  card  .  56 

12.  The  OPTl  card .  57 

13.  The  SCN3  card .  57 

14.  The  TCAL  card .  58 


5 


15.  TlieSVARcard  .  58 

16.  The  GO  card  .  59 

17.  The  DONE  card  .  60 

18.  Run  times  for  selected  run  conditions .  60 

Appendix  Table 

C-1.  Limitations  caused  by  model  assumptions . 100 


6 


1.  Introduction 


The  Fire-Induced  Transmission  and  Turbulence  Effects  (FITTE)  module  predicts 
transmittance  through  fire  plumes,  emitted  radiance  fi-om  fires  and  fire  plumes, 
and,  optionally,  effects  of  fire  and  fire  plume  turbulence  on  laser  propagation 
for  a  given  line  of  sight  (LOS).  The  FGLOW  model  option  extends  the 
calculations  to  multiple  LOSs  to  predict  the  radiant  image  of  a  fire  or  fire 
plume  segment  seen  by  an  imaging  system.  The  fires  represent  localized 
sources  of  burning  diesel  fuel,  motor  oil,  and  rubber  (simulated  burning 
vehicles).  FITTE  also  predicts  the  path-integrated  particulate  concentration, 
and,  if  a  target  temperature  is  specified,  the  attenuated  thermal  radiance  from 
the  target  at  the  observer  position.  If  the  calculation  is  performed  for  a  single 
wavelength,  the  model  predicts  beam  spread  and  wander  for  a  laser  beam  of 
that  wavelength. 

The  FGLOW  option  performs  calculations  for  a  fan-shaped  rectangular  set  of 
LOSs.  It  creates  a  file  of  emitted  radiance  values,  that  represent  the  radiant 
image  of  the  fire  or  fire  plume  seen  by  an  imaging  system,  and  of  transmittance 
values  that  represent  the  attenuation  of  file  background  radiance.  An  option 
allows  output  of  an  apparent  temperature  image  of  the  fire  or  fire  plume. 

The  information  in  the  technical  manual  (sections  1  through  3)  describes  the 
model  at  the  time  it  was  furnished  to  the  Electro-Optical  Systems  Atmospheric 
Effects  Library  (EOSAEL).  Changes  at  future  dates  will  be  documented  in  a 
README  or  similar  text  file,  which  should  be  distributed  with  the  code. 
Similar  comments  apply  to  the  user’s  manual  (sections  4  and  5). 


1.1  Names  of  Things 

EOSAEL  standard  typefaces  are  used  for  descriptions: 

1.  Names  of  modules  are  in  univers  style. 

2.  Variable  and  subroutine  names  are  in  courier  style. 

3.  Sample  input  and  output  are  in  typewriter  style,  monospaced  to  allow 
column  alignment. 
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1.2  Availability 

EOSAEL  92  is  available  at  no  cost  to  U.S.  Government  agencies,  specified 
Allied  organizations,  and  their  authorized  contractors.  U.S.  Government 
agencies  needing  EOSAEL  92  should  send  a  letter  of  request,  signed  by  a 
branch  chief  or  division  director,  to  U.S.  Army  Research  Laboratory  (ARL). 
Contractors  should  have  their  Government  contract  monitor  send  the  letter  of 
request.  Allied  organizations  must  request  EOSAEL  92  through  their  national 
representative.  The  EOSAEL  92  point  of  contact  at  ARL  is  Dr.  Alan  Wetmore. 

Include  intended  uses  with  requests.  Also,  include  the  type  of  nine-track  tape 
necessary  for  computer  execution.  Tape  formats  are  as  follows;  ASCII,  UNIX 
tar  format  in  1600  or  6250  bpi,  or  SUN  cartridge.  EOSAEL  92  cannot  be 
supplied  on  other  media.  Documentation  for  modules  is  included. 

1.2.1  Mailing  Address 
Director 

Battlefield  Environment  Directorate 
Army  Research  Laboratory 
AMSRL-BE-S  (Attn:  Dr.  A.  Wetmore) 

White  Sands  Missile  Range,  NM  88002-5501 

1.2.2  Phone  and  Electronic  Mail 

(505)  678-5563 
FAX  (505)  678-2432 

DSN  258-5563 

email  awetmore@arl.mil 
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2.  Background 

The  FITTE/FGLOW  model  was  developed  to  provide  relatively  rapid 
calculations  of  the  effects  of  battlefield  fires  on  detection,  laser-designator,  and 
weapon  systems.  Earlier  versions  of  the  models  treated  the  fire  and  fire  plume 
as  static  mean-value  entities. 

Although  much  work  has  been  done  on  modeling  laser  propagation  in  the 
atmosphere,  the  situations  modeled  were  those  for  which  homogeneous, 
isotropic  turbulence  could  be  assumed.  The  assumptions  permit  simplifications 
in  calculations  of  short-  and  long-time  average  spot  size  and  average  centroid 
jitter  not  valid  for  propagation  through  a  fire  plume. 

Fire  plumes  are  regions  of  space  where  there  are  rapidly  changing  atmospheric 
conditions.  Special  algorithms  are  needed  to  treat  the  effects  of  fire  plume 
turbulence  on  propagation  because  of  the  rapid  changes  in  refi'active  index 
within  the  plume,  the  large  and  nonuniform  increase  in  the  refractive  index 
structure  constant,  and  the  fact  that  the  turbulence  effects  depend  on  path- 
weighted  turbulence  parameters. 

A  preliminary  model  of  the  effects  of  a  vertical  fire  plume  on  laser  propagation 
was  developed  (Thompson  and  DeVore  1981)  for  the  U.S.  Army  Atmospheric 
Sciences  Laboratory  (ASL).  This  model  was  modified  (Bruce,  Sutherland,  and 
Larson  1982)  to  handle  bent-over  plumes  with  obscuration  fi-om  smoke  for 
EOSAEL  82.  The  model  predicted  transmittance  throu^  the  plume  for  a  given 
LOS  and  wavelength,  and  it  optionally  predicted  magnitudes  of  beam  wander 
and  beam  spread.  The  transmittance  depended  only  on  the  plume  aerosols 
because  the  major  gaseous  effluents,  water  vapor  and  carbon  dioxide,  are 
responsible  for  less  than  5  percent  of  the  total  extinction  in  the  visual  and 
infi'ared  (IR)  atmospheric  window  regions  of  the  spectrum  where  many  Army 
systems  operate. 

Fire  and  plume  parameters  (fire  temperature  and  plume  vertical  velocity)  were 
obtained  fi-om  laboratory  and  field  measurements  (D.  Bruce,  unpublished  data; 
C.  W.  Bruce,  private  communication  1982).  Smoke  optical  parameters  were 
obtained  by  a  semiempirical  approach.  A  theoretical  approach  to  determining 
specific  extinction  coefficients  for  carbonaceous  aerosols  has  been  described 
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(Chylek  et  al.  1981).  This  entails  assuming  equivalent  spheres  with  appropriate 
particle-to-void  ratios  and  applying  Mie  scattering  theory.  For  the  FITTE 
model,  particle-to-void  theory  was  modified  to  give  agreement  with 
transmissometer  data  (Sutherland  and  Walker  1982;  Sutherland,  Hoock,  and 
Khanna  1984). 

Minor  changes  were  made  in  the  FITTE  model  for  EOSAEL  84.  The  changes, 
primarily,  were  to  the  calculation  of  radiative  transfer.  Model  output  was 
expanded  to  include  attenuated  target  radiance  from  a  hot  target  and  emitted 
radiance  from  the  plume,  which  made  it  possible  to  predict  target  contrast. 

The  Forschungsinstitut  fiir  Optik  commissioned  the  development  of  the  FEU  ER 
code  (Manning  1985),  which  was  a  modification  of  the  1984  version  of  FITTE 
that  added  effects  of  absorption  and  emission  by  the  major  effluent  gases  to  the 
predictions  and  the  capability  to  perform  calculations  over  a  waveband.  The 
work  was  based  on  a  modification  of  code  from  the  ATLES  model  (Young 
1977).  Temperature-dependent  values  of  absorption  coefficients  are  included 
for  water  vapor  in  the  spectral  region  from  2.222  to  200.000  pm  (4500  to 
50  cm  '),  and  for  carbon  dioxide  in  the  regions  from  2.702  to  3.225  pm  (3700 
to  3100  cm"')  and  4.160  to  5.000  pm  (2400  to  2000  cm'').  Data  are  included 
for  temperatures  from  100  to  3000  K.  The  FEUER  code  was  supplied  to  ASL 
for  incorporation  into  the  EOSAEL  87  version  of  FITTE. 

The  FGLOW  model  (Bruce  1988)  also  was  added  for  EOSAEL  87.  FGLOW 
is  an  extension  of  FITTE  that  performs  transmittance  and  emitted  radiance 
calculations  for  a  set  of  LOSs  in  a  fan-shaped  rectangular  array.  The  output 
can  be  used  to  predict  the  mean  image  of  the  fire  or  fire  plume  that  would  be 
seen  by  an  imaging  sensor,  or  it  may  be  used  in  conjunction  with  an  appropriate 
algorithm  to  add  a  fire  plume  to  a  background  image. 

Significant  changes  were  made  for  the  first  revision  of  EOSAEL  87  FITTE. 
Some  changes  result  from  work  done  (to  increase  the  speed  of  calculation)  for 
the  Forschungsinstitut  fur  Optik  by  OptiMetrics,  Inc.  (Manning  and 
Kebschull  1988),  and  other  changes  result  from  work  at  ARL.  Most  of  the 
changes  are  transparent  to  the  user.  They  cause  an  increase  in  the  speed  of 
calculations  by  up  to  a  factor  of  10,  particularly  for  waveband  and  FGLOW 
image  calculations. 


Time-dependent  spatial  variations  were  added  to  the  FITTE/FGLOW  model  for 
EOSAEL  92.  The  variations  were  accomplished  by  adding  time-dependent 
fluctuations  to  the  centerline  position  and  to  most  of  the  mean-value  Gaussian- 
distributed  parameters.  Other  changes  include  options  to  perform  time-series 
calculations  and  to  vary  one  or  more  of  the  fire  parameters. 

2.1  Conceptual  Overview  of  FITTE 

Simplified  physical  models  of  fires  and  fire  plumes  are  combined  with  a  model 
for  radiative  transfer  to  predict  the  effects  of  battlefield  fire  plumes  on 
propagation  of  visual  and  IR  radiation  or  to  predict  values  needed  to  evaluate 
the  effects  of  fire  plumes  on  imaging  sensors.  A  specially  tailored  model  of  the 
turbulence  structure  function  within  a  fire  plume  is  used  to  predict  beam  spread 
and  wander  for  laser  propagation.  Empirical  data  and  results  of  theoretical 
calculations  are  used  for  fire  and  aerosol  parameters. 

2.2  The  Physical  Models 

2.2.1  The  Fire  and  Fire  Plume  Models 

The  fire  and  fire  plume  are  described  in  terms  of  mean  time-averaged 
parameters  with  time-varying  fluctuations  superimposed  to  create  dynamical 
variation.  The  description  of  the  mean  value  model  is  given  first.  A 
description  of  the  fluctuation  modeling  follows. 

2. 2. 1. 1  The  Time-Averaged  Fire. — FITTE/FGLOW  uses  an  extremely  simple  model  for 

the  fire  itself  A  fire  is  modeled  as  a  circular  disk  with  fixed  mean  temperature 
r„,  radius  R^,  and  bum  time  4uni»  and  specific  quantities  of  fuel.  Observational 
evidence  shows  that  the  bum  time  should  be  a  function  of  windspeed,  and  such 
dependence  may  be  incorporated  into  the  model  after  quantitative  data  become 
available. 

The  fire  is  centered  at  the  origin  of  the  model  coordinate  system.  The  model 
coordinate  system  is  a  right-handed  coordinate  system  with  the  positive  x-axis 
coinciding  with  the  downwind  direction.  The  coordinate  system  is  used  because 
it  affords  the  simplest  possible  equations.  Transformation  to  and  fi-om  other 
coordinate  systems  is  done  using  standard  EOSAEL  conventions. 
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The  user  specifies  a  fire  type  firom  several  types  that  represent  burning  vehicles 
or  vehicular  fuels.  The  model  is  restricted  to  these  fires  because  they  are  the 
only  type  of  fiiel  mixture  for  which  a  relatively  complete  set  of  data  is 
available.  The  fire  size  is  restricted  so  the  fire  does  not  significantly  change  the 
atmospheric  parameters  used  in  the  calculations. 

The  fire  characteristics  were  chosen  based  on  data  fi-om  Battlefield  Induced 
Contamination  Tests  I  and  III  (BICT  I  and  BICT  III)  (Kennedy  1981;  Kennedy 
1982).  The  characteristics  were  designed  to  give  realistic  values  for  mean  heat 
and  aerosol  fluxes.  Types  1,  2,  and  3  correspond  roughly  to  small 
medium  (truck),  and  large  (tank)  burning  vehicles,  while  type  4  is  a  larger  fire 
that  represents  a  small  fuel  depot.  The  four  types  should  be  reasonable 
representations  of  battlefield  fires.  Table  1  lists  defining  values  for  the  fires. 


Table  1.  Values  defining  model  fires 


Type 

Diesel  Fuel 

Motor  Oil 

Rubber 

Radius 

(gal) 

[kg] 

(qt) 

[kg] 

(kg) 

(m) 

Jeep 

1 

5 

[16.821] 

6 

[5.0463] 

40 

0.6 

Truck 

2 

20 

[67.284] 

8 

[6.7284] 

100 

1.0 

Tank 

3 

125 

[420.525] 

16 

[13.4568] 

250 

2.0 

Fuel  Depot 

4 

800 

[2,691.360] 

160 

[134.5680] 

0 

4.0 

Quantities  derived  from  the  source  parameter  values  are  the  aerosol  (smoke) 
mass  emission  rate,  the  heat  emission  rate,  and  the  number  densities  of  the 
major  molecular  combustion  products  (carbon  dioxide  and  water  vapor).  The 
values  are  converted  to  fluxes  through  the  plume  cross-sectional  area. 

Calculation  of  the  aerosol  mass  emission  rate  Qy^  requires  a  value  for  the 
particulate  emission  factor  the  ratio  of  the  mass  of  particulates  produced  to 
the  mass  of  fuel  consumed  in  the  combustion  process.  This  quantity  is  difficult 
to  determine  and  is  estimated  fi-om  experimental  data  (Sutherland  and  Walker 
1982)  to  be  0.160  g/g.  The  mass  emission  rate  is  modeled  as 
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(1) 


Mdf  +  ^RU 


''bum 


where  M  is  the  total  mass  (in  grams)  of  diesel  fuel  {DF),  motor  oil  {MO),  or 
rubber  {RU),  and  tbum  is  the  total  bum  time,  set  to  25  min  to  obtain  heat  and 
mass  flux  levels  similar  to  those  at  the  BICT  tests. 

The  heat  emission  rate  is  given  by 

Q  =(!-£)  ^  ^  (2) 

^bum 

where  the  .E’s  in  the  numerator  are  the  heat  of  combustion  of  the  various  fuels. 

The  factor  (1  -  arises  to  account  for  the  fact  that  particulates  are  unbumed 
residue  of  the  combustion  process.  Values  for  the  heats  of  combustion  are 
summarized  (Gomez  et  al.  1980)  from  data  (Turner,  Eitner,  and  Manning  1980; 
Levitt  and  Levitt  1982). 

The  molecular  concentrations  (number  of  molecules/cc)  of  the  major 
combustion  products  are  arrived  at  as  follows.  Combustion  of  diesel  fuel  is 
represented  (Leslie  et  al.  1980)  by 

(3) 

2  g  +  2.9  Oj  +  10.93  -  2  CO^  +  1.8  Hp  +  10.93  N^. 


The  reaction  in  equation  (3)  is  a  composite  representation  of  a  set  of  reactions, 
one  for  each  of  the  various  fuel  molecules.  Although  the  reaction  in 
equation  (3),  generally,  is  used  to  approximate  combustion  in  an  engine,  it  is 
equally  valid  for  combustion  in  air.  The  molecular  concentration  of  carbon 
dioxide  caused  by  combustion,  Csj,  in  a  1-cm  high  volume  above  a  fire  of  area 
ap  cm^  is 
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1 


1 


'SI 


hydrocarbon  fuel  weight 
molecular  weight  of  CH.^  g 


(4) 


where  Nj^is  Avogadro’s  number  and  is  the  effluent  velocity  (cm/s)  just  above 
the  fire. 

From  equation  (3),  the  molecular  concentration  of  water  vapor,  0^2,  produced 
by  combustion  is  0.9  times  that  of  carbon  dioxide. 

The  mean  temperature  (923.16  K)  for  the  top  hat  fire  distribution  was  chosen 
based  on  unpublished  data  taken  by  the  author  during  BICT  I  (Kennedy  1981). 
This  translates  to  a  peak  temperature  of  roughly  2000  K  for  the  equivalent 
Gaussian  distribution  used  in  the  propagation  model. 

2.2. 1.2  The  Time-Averaged  Fire  Plume. — FITTE  uses  a  top  hat  plume  model  to 
describe  the  transport  of  general  curved  plumes  as  well  as  the  special  cases  of 
horizontal  and  vertical  plumes.  The  top  hat  distribution  has  a  single  value 
within  the  plume  radius  and  is  zero  outside  the  plume. 

The  two-thirds  law  model  (Briggs  1 969)  is  used  to  predict  the  time-averaged 
plume  centerline.  The  plume  centerline  equation  for  a  finite  source  at  the 
origin  of  the  coordinate  system  is 

z  +  Zy  =  ^  ix  +  (5) 


where  x  and  z  are  plume  centerline  coordinates  in  the  x-z  plane,  Xy  and  Zy  are 
constant  coordinates  of  a  virtual  point  source  and  are  chosen  to  match  boundary 
conditions  at  the  fire  source,  and  U  is  the  ambient  mean  windspeed  (assumed 
uniform  with  height). 

U  is  not  permitted  to  be  less  than  0.10  m/s  to  ensure  that  equation  (5)  is  well 
defined.  Figure  1  shows  a  time-averaged  plume  in  the  model  coordinate 
system. 
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Figure  1.  Side  view  contour  plot  of  a  fire  plume  in  the  model  coordinate  system. 
The  plume  edges  shown  are  at  the  at  the  e'^  points  of  a  conical  Gaussian  distribution. 
The  centerline  is  sketched. 


The  constant  C  is  given  by  (Weil  1981)  as 


where  ^  is  an  entrainment  constant;  a  constant  that  describes  the  fractional 
increase  in  plume  volume  per  unit  distance  along  the  centerline  caused  by 
diffrision,  and  F  is  the  buoyancy  flux. 

F  is  defined  as 

F  =  ~  (7) 
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where  g  is  the  gravitational  acceleration,  is  the  source  heat  flux,  Cp  is  the 
specific  heat  of  air  at  constant  pressure,  is  the  ambient  air  density,  is  the 
ambient  temperature,  is  the  source  temperature,  is  the  initial  plume 
vertical  velocity,  and  is  the  source  radius. 

If  it  is  assumed  that  the  buoyancy  flux  is  conserved,  equation  (7)  applies  at  any 
height  z  if  the  values  and  are  replaced  by  T(z),  W(z),  and  R(z), 

respectively. 


Using  a  geometrically  based  argument  similar  to  those  in  the  literature  (Hoult, 
Fay,  and  Forney  1969),  the  entrainment  parameter,  for  a  general  plume  can 
be  expressed  as 


^  =  n 


1  - 


V  / 


+  u 


u 

\  ^  I 


(8) 


where  tj  is  the  empirically  determined  entrainment  constant  for  a  vertical  plume, 
V  is  the  empirically  determined  entrainment  constant  for  a  horizontal  plume,  and 
S  is  the  plume  velocity  along  the  centerline. 

Experimentally  determined  values  used  for  entrainment  into  vertical  and 
horizontal  plumes  are  r;  =  0.12  and  v  =  0.60,  respectively  (Hoult  and 
Weil  1972),  although  values  have  been  reported  (Briggs  1969)  as  low  as  0.07 
for  vertical  plumes  and  as  high  as  1.0  for  horizontal  plumes.  ^  must  be 
constant  to  maintain  model  simplicity,  so  the  value  of  plume  velocity  at  the 
origin  S(0)  is  used  in  the  defining  equation. 

The  mean  velocity  along  the  centerline  S  is  assumed  to  be  the  vector  sum  of  the 
ambient  mean  wind  vector  U  and  the  mean  vertical  plume  velocity  W.  S  and 
W  may  be  written  as  functions  of  z  or,  alternatively,  in  terms  of  D,  the  distance 
from  the  origin  along  the  centerline.  The  time-averaged  plume  parameters, 
assumed  to  have  Gaussian  distributions  about  the  plume  centerline  (George, 
Alpert,  and  Tamanini  1977),  are  most  easily  expressed  as  functions  of  D  and 
of  the  plume  width  parameter  g(D). 

The  relation  between  the  value  of  the  centerline  coordinate  z  and  the  distance 
D  of  that  point  from  the  origin  is  obtained  as  follows.  The  ambient  wind  is 
assumed  to  be  horizontal,  so  W(z)  can  be  expressed  as 


(9) 


W(z)  ^  jdzldt)  ^  & 
U  (dx/dt)  dx 


at  any  point  on  the  centerline. 


Substitution  for  dz/dx  from  the  plume  centerline  equation  gives 


(10) 


Because  the  vertical  wind  velocity  at  the  origin  is  known  and  ^  is  defined 
in  terms  of  known  quantities,  the  constants  Zy,  Xy,  and  C  can  be  found. 


The  plume  model  assumes  that  the  time-averaged  plume  parameter  distributions 
vary  as  functions  of  D.  Therefore,  it  is  usually  more  convenient  to  express  the 
parameters  as  functions  of  D  rather  than  of  centerline  height  z.  The  relationship 
between  D  and  z  is  found  by  integrating 


dD  =  )j  {dx/dz)^  +  1  dz 


from  the  origin  to  the  point  of  interest,  z^i,  along  the  centerline. 
Integrating  equation  (11)  gives 


2.2. 1.3  Time-Averaged  Distribution  Functions. — Spatially  resolved  measurements  of 
time-averaged  plume  parameters,  such  as  temperature  and  vertical  velocity,  can 
be  described  by  circularly  symmetrical  Gaussian  distributions  at  fixed  distances 
from  the  fire  (Hoult,  Fay,  and  Forney  1969).  Such  distributions  are  used  to 
describe  the  mean  values  of  the  temperature,  effluent  concentrations  (of  aerosols 
and  molecular  combustion  products),  and  the  vertical  velocity  of  the  plume 
when  calculations  are  performed  to  predict  the  transmittance  and  turbulence 
effects.  For  a  vertical  plume,  such  a  distribution  has  the  form 
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(13) 


A^{r,z)  =  A^{Q,z)'e 

where  A^(...)  is  the  parameter  amplitude,  r  is  the  radial  distance  from  the 
centerline,  and  ct(z)  is  the  diffrision  function  as  a  flmction  of  distance  from  the 
fire. 

The  distribution  function  behavior  remains  the  same  for  the  general  curved 
plume  used  in  FITTE,  but  the  mathematical  description  becomes  more  complex. 
Conservation  assumptions  determine  the  decrease  in  peak  amplitude  of  the 
parameters  as  a  function  of  distance  along  the  centerline.  The  general  variation 
of  the  plume  distribution  flmctions  is  shown  in  figure  2. 
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The  difiusion  function  a  is  related  to  the  entrainment  parameter  by 

(14) 


a(D)  =  o(0)  + 


6.06 


and  to  the  plume  radius  R  by 


R  =  2.0  a(D). 


(15) 


The  definition  of  the  plume  radius  agrees  rougjily  with  the  apparent  plume  edge 
seen  by  an  observer.  In  performing  the  propagation  integrals,  FITTE  uses  the 
plume  parameter  distribution  functions  when  a  LOS  point  is  within  3a  of  the 
centerline  (when  the  parameter  is  greater  than  approximately  1  percent  of  its 
peak  value).  The  same  radius  is  used  for  all  the  plume  parameters. 

The  relation  between  equivalent  top  hat  and  Gaussian  plume  distributions  has 
been  described  elsewhere  (Thompson  and  DeVore  1981).  For  equivalent 
distributions,  the  integrated  top  hat  values  of  mass,  momentum,  and  buoyancy 
are  required  to  equal  truncated  integrals  of  the  corresponding  Gaussian  plume 
parameters.  Figure  3  shows  equivalent  distributions,  with  the  shaded  portions 
in  the  left  half  indicating  the  areas  mapped  by  the  transformation  between  the 
distributions.  The  maximum  of  the  Gaussian  distribution  is  approximately  2.3 1 
times  the  top  hat  amplitude  for  the  effluent  densities.  The  relation  for  the 
excess  temperature  above  ambient  is  more  complex,  and  the  ratio  ranges  from 
about  2.3  to  2.9. 


Figure  3.  Equivalent  top  hat  and  Gaussian 
distributions. 
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The  equation  used  to  describe  effluent  aerosol  concentration  is  based  on  the 
equation  for  a  classic  horizontal  Gaussian  plume; 


X(  x,y,z) 


g  2a^ 


(z  -  Zc)  _  *  ^c) 

2a2  +  R^e 


(16) 


where  x  is  the  aerosol  concentration  at  any  point  specified  in  the  model 
coordinate  system,  Qy^  is  the  aerosol  mass  emission  rate,  U  is  the  mean  ambient 
wind  velocity,  o(x)  is  the  diffusion  function,  and  Zq  is  the  height  of  the  plume 
centerline  at  x. 

The  effluent  molecular  combustion  product  concentrations  are  described  in  a 
similar  fashion. 


Equation  (16)  applies  only  in  the  extreme  case  of  a  horizontal  plume.  To 
include  plume  curvature  or  a  vertical  plume,  the  equation  must  be  modified  so 
the  parameters  are  related  to  the  direction  of  the  mean  flow  at  any  point  along 
the  centerline.  This  transformation,  although  simple  in  principle,  introduces 
geometrical  and  mathematical  complexities,  because  the  direction  of  the  mean 
flow  vector  is  a  function  of  location  for  a  curved  plume. 

For  the  general  curved  plume,  the  analog  of  equation  (16)  is 


x(A  K  y) 


2-31 
4tz  Sa^ 


K^e  +  K^e 


(  2hc  +  K^h  f 
2? 


(17) 


where  D  is  the  distance  along  the  plume  centerline  to  a  point  (xc,,  0,  z^),  h  is 
the  distance  (measured  perpendicular  to  the  mean  flow  vector  S  at  D)  from  the 
plume  centerline  to  the  projection  of  a  point  (xl,  0,  yi)  on  the  LOS  onto  the  x-z 
plane,  and  is  the  distance  of  the  plume  centerline  from  the  ground  (also 
measured  perpendicular  to  S  at  D). 

ct(D)  is  related  to  the  plume  radius  at  D,  as  illustrated  in  figure  2.  The 
coordinate  y  is  the  distance  from  the  plume  centerline  perpendicular  to  the  x-z 
plane. 


The  constants  K^,  K2,  and  in  equation  (17)  take  on  values  appropriate  for 
specific  situations.  For  a  horizontal  plume,  =  K2  =  1.0.  For  a  vertical 
plume,  K^  =  I  and  ^^2  =  0.  The  constants  account  for  the  reflection  from  the 
ground  for  quasihorizontal  plumes.  The  conditions  under  which  the  model 
changes  from  one  set  of  values  to  the  other  for  bent-over  plumes  is  described 
later.  To  obtain  the  appropriate  value  for  the  reflected  concentration  when  there 
is  reflection,  is  set  equal  to  +1  if  the  point  at  which  the  function  is  to  be 
evaluated  lies  above  the  centerline,  or  to  -1  if  the  point  lies  below  the 
centerline. 

An  accurate  description  of  the  plume  temperature  distribution  is  needed  for 
calculations  of  the  turbulence  effects.  If  values  from  the  equation  for  pressure 
equilibrium. 


p(r)-f(r)  = 

are  substituted  into  the  equation  for  the  equivalence  of  mass, 

Rtu 

=  j  2%p{r)  rdr, 

0 


(18) 


(19) 


the  centerline  temperature  increment  above  ambient  can  be  determined  for  an 
equivalent  truncated  Gaussian  plume  (one  for  which  the  integral  of  the  top  hat 
distribution  is  equal  to  the  integral  of  the  truncated  Gaussian  distribution). 
Overbars  indicate  the  mean  values  associated  with  the  top  hat  model. 
Experiments  indicate  that  equivalence  of  the  distributions  should  be  defined  for 
a  top  hat  radius  i?TH  =  2.0  a(D),  which  corresponds  roughly  to  the  visible  edge 
of  the  plume.  This  definition  leads  to  the  expression  for  the  centerline 
temperature: 


Tciz)  = 


(20) 
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The  relations  between  the  peak  and  mean  values  of  the  effluent  concentration 
and  the  vertical  velocity  are  found  in  a  similar  manner.  It  is  assumed  that  all 
of  these  distributions  are  characterized  by  the  same  plume  radius.  Experiments 
(D.  Bruce,  unpublished  data)  indicate  that  this  assumption  is  approximately 
correct. 


Finally,  the  gradient  of  the  amplitude  of  the  plume  temperature  distribution 
must  be  found  as  a  function  of  distance  D  along  the  centerline  (or,  alternatively, 
as  functions  of  the  centerline  height  z)  for  use  in  the  turbulence  effects 
calculations.  This  is  done  by  expressing  the  centerline  temperature  as  a 
function  of  z  and  taking  the  derivative  with  respect  to  z  to  obtain 


dTg  _  dT 

dz  df  dz 


where 

df 


2 


T 


IT.  /  T 

e 


2T,  I  T 

e 


(21a) 


(21b) 


^  ^  y,  p(  2C^>^  Z~^'^  +  eiU^'^WR  ^ 

dz  {F  -  gWR'^  j 


(21c) 


R  and  W  are  treated  as  explicit  functions  of  z  in  arriving  at  equation  (21c)  and 
are  treated  as  such  in  its  evaluation. 

The  change  in  amplitude  of  the  effluent  flux  distributions  as  functions  of  D  are 
found  by  observing  that  the  mass  flow  through  all  plume  cross  sections  must 
remain  constant  to  prevent  buildup  of  mass  at  any  location,  so  at  any  distance 
D  from  the  origin,  the  effluent  flux  must  be  equal  to  the  initial  flux  times  the 
product  of  the  plume  velocity  S  and  the  cross-sectional  area  at  the  origin 
divided  by  the  same  product  evaluated  at  distance  D. 

2. 2. 1.4  Aerosol  Optical  Properties. — ^Visual  and  IR  aerosol  optical  parameters  were 

based  on  theoretical  calculations  modified  to  give  agreement  with 
transmissometer  data  from  BICT  I.  The  aerosol-specific  extinction  coefficients 
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a  for  various  wavelengths  were  based  on  Mie  scattering  theory  calculations. 
The  calculations  used  a  theory  (Chylek  et  al.  1981)  that  assumes  carbon-based 
aerosols  can  be  replaced  by  equivalent  spheres  with  appropriate  particle-to-void 
ratios.  The  theory  was  modified  to  give  agreement  with  ratios  of  measured 
transmittance  for  different  wavelengths.  Details  are  given  elsewhere 
(Sutherland  and  Walker  1982;  Sutherland,  Hoock,  and  Khana  1984)  and  will 
not  be  repeated  here. 

The  values  of  single  scattering  albedo  to  were  derived  from  measurements  of 
specific  extinction  and  absorption  (C.  W.  Bruce,  private  communication  1981) 
or  were  taken  from  Mie  calculations  (Khanna  and  Ammon  1983)  for 
aerosol-size  distributions  and  porosity  characteristics  based  on  measurements 
(Pinnick  et  al.  1982). 

Table  2  lists  the  values  used  in  FITTE.  Interpolation  is  used  to  calculate  the 
specific  extinction  coefficient  and  albedo  at  intermediate  wavelengths.  The 
10.6-|j,m  and  millimeter  wavelength  values  are  from  measurements  (C.  W. 
Bruce,  private  communication  1982). 


Table  2.  Particulate  optical  parameters 


Wavelength 

(pm) 

Specific  Extinction 
Coefficient 
(mVg) 

Single  Scattering 

Albedo 

0.40 

7.41 

0.25150 

0.55 

5.57 

0.25000 

1.06 

2.94 

0.24493 

1.50 

2.23 

0.24055 

2.50 

1.61 

0.23060 

3.80 

1.34 

0.21766 

6.00 

1.14 

0.19577 

8.50 

1.01 

0.17090 

10.60 

1.00 

0.15000 

8570.00 

0.002 

0.00000 
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2.2. 1.5  Fire  Plume  Fluctuations. — ^The  space-time  variation  of  the  fire  plume  is 
handled  as  a  two-part  model.  The  fluctuations  of  the  plume  about  the  mean 
centerline  are  treated  as  functions  of  windspeed  and  wind-direction  fluctuations, 
while  the  fluctuations  of  the  plume  parameters  are  based  on  the  assumption  that 
they  are  dependent  on  turbulent  mixing. 

A  time-history  data  file  FIT4D.DAT  is  used  to  store  the  information  needed  for 
the  fluctuation  model  in  FITTE/FGLOW.  The  FIT4D.DAT  file  contains 
2000  rows  and  3  columns  of  data.  The  rows  of  data  have  a  1-s  time  increment 
and  represent  a  time  history  of  wind  fluctuations  since  the  start  of  the  fire. 
Column  1  contains  a  random  number  (range:  -0.5  to  0.5)  used  as  input  for 
calculating  plume  parameter  fluctuations  at  a  given  time,  column  2  contains  the 
normalized  cumulative  fluctuation  of  the  downwind  (x-)  component  of  the 
wind,  and  column  3  contains  the  normalized  cumulative  fluctuation  of  the 
crosswind  (y-)  component  of  the  wind.  The  model  was  developed  and  tested 
with  simulated  data.  Files  developed  from  actual  wind  data  as  described  in 
appendix  A  can  be  used  after  model  operation  is  checked.  The  file  structure  is 
discussed  further  in  appendix  A,  which  gives  information  for  creating  a  time- 
history  file  from  wind  data. 

When  FITTE  is  run,  a  current  time  is  assigned.  The  current  time  may  be  the 
default  time  (100  s  after  the  start  of  the  fire)  or  another  time  read  from  an  input 
file. 

The  displacement  of  any  point  P  on  the  centerline  from  its  mean  position  at  the 
current  time  is  calculated  as  follows:  The  distance  along  the  mean  centerline 
from  the  origin  to  point  P  is  calculated  using  equation  (12).  The  average  plume 
velocity  S  to  reach  point  P  is  half  the  sum  of  S(P)  and  S(0).  The  ratio  of 
distance  to  mean  velocity  gives  the  time  taken  to  reach  point  P  from  the  fire. 
The  time  to  reach  point  P  from  the  fire  is  used  to  assign  a  source  time  (the  time 
when  the  plume  segment  currently  at  point  P  left  the  source).  The  current  time 
and  source  time  are  used  as  pointers  to  rows  of  the  data  array.  Using  the 
second  (or  third)  column  of  the  array,  the  difference  in  value  between  the 
values  in  these  rows  multiplied  by  the  mean  windspeed  U  gives  the  x-  (  or  y-) 
displacement  of  point  P  from  the  mean  centerline  position.  This  permits 
calculation  of  the  current  position  of  point  P  to  determine  the  plume  parameter 
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distributions  at  points  on  an  LOS  near  point  P.  Interpolation  takes  place  when 
the  time  to  reach  point  P  is  not  an  integer. 

The  random  numbers  in  the  first  column  of  the  FIT4D.DAT  file  are  used  in  two 
ways: 

(1)  Assign  scaled  firactional  fluctuations  to  the  plume  radius  and  to  the  heat  and 
effluent  production.  The  scaling  limits  the  fluctuations  to  no  more  than 
5  percent  of  the  mean  value  of  the  corresponding  parameter. 

(2)  The  numbers  are  used  as  input  for  calculation  of  the  fluctuations  of  the 
aerosol  and  effluent  gas  concentrations  and  of  the  plume  temperature. 

Modeling  of  the  plume  parameter  fluctuations  is  biised  on  the  assumption  that 
they  are  driven  by  atmospheric  turbulence.  Therefore,  they  should  exhibit  the 
same  power  spectral  density  distribution  as  found  for  atmospheric  parameters 
within  an  appropriate  fi’equency  range.  This  is  achieved  by  the  choice  of  the 
coefficients  used  for  the  series  that  describes  the  fluctuations.  The  fluctuations 
are  represented  by  a  series  with  the  following  form: 


6  =  ^  cos(  b^x  +  (pjj. )  . 

1 

The  phase  angles  (p^  for  the  terms  are  obtained  fi'om  the  random  numbers  stored 
in  the  first  column  of  the  FIT4D.DAT  file.  Figure  4  is  an  example  of  a  space- 
and  time-varying  (4D)  plume. 

2.2,2  The  Radiative  Transfer  Model 

The  radiative  transfer  equation  for  a  medium  in  which  aerosol  scattering  and 
absorption  and  molecular  absorption  occur  is 

^  -  G/  +  (K^  +  (23) 
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show  5-m  spacing.  The  dashed  line  is  the  sketched  centerline. 


where  /  is  the  local  spectral  radiance  (W/mVsr/cm),  s  is  the  distance  along  the 
LOS,  e  is  the  total  extinction  coefficient  (m'*),  #c^  and  are  the  aerosol  and 
molecular  absorption  coefficients  (m'‘),  and  B  is  the  blackbody  function. 


BiX,  T)  = 


Inhc^ 


1 


ghclkkT  _  2 


(24) 


where  h  is  Planck’s  constant,  k  is  Boltzmann’s  constant,  and  c  is  the  speed  of 
li^t  in  vacuum. 

The  radiance  /  has  units  of  W/mVsr/)im,  and  the  blackbody  irradiance  B  has 
units  of  W/mV(im. 

For  the  aerosol,  the  absorption  coefficient  is  related  to  the  extinction 
coefficient  e  by 
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=  (1  -  o))e  =  (1  -  co)ax 


(25) 


where  o:  is  the  specific  mass  extinction  coefficient  of  the  aerosol  (mVg),  oi  is  the 
single-scattering  albedo,  and  x  is  the  local  aerosol  concentration  (g/m^). 

In  FITTE,  only  the  effluent  aerosol  and  the  major  molecular  combustion 
products  (water  vapor  and  carbon  dioxide)  are  considered  in  the  calculation. 
The  aerosols  are  the  dominant  species  for  calculations  in  the  visible  and  IR 
atmospheric  window  regions  (Sutherland  and  Walker  1982;  IChanna,  Burlbaw, 
and  Deepak  1982).  If  the  path  from  target  to  observer  is  treated  as  a  series  of 
segments  for  which  the  effluent  concentrations  and  temperature  are  nearly 
homogeneous,  the  radiative  transfer  equation  can  be  integrated  for  each  path 
segment  to  give 


I  =  Le 


ne  '■ 


(1  -  o))ax  +  5(1  -  e~^)  ] 


(26) 


where  r  is  the  optical  depth  of  the  segment,  I  is  the  spectral  radiance  in  the 
segment,  and  /g  is  the  spectral  radiance  entering  the  segment  from  all  previous 
segments. 

The  radiative  transfer  calculation  is  performed  by  evaluating  this  expression  for 
each  segment  from  the  target  to  the  observer  using  /  for  each  segment  and  /g 
for  the  following  segment. 

The  extinction  coefficients  are  functions  of  wavelength  and  temperature.  The 
values  used  for  the  aerosol  are  temperature-average  values  from  measurements; 
the  values  for  the  molecular  combustion  products  are  temperature  dependent  and 
are  taken  from  the  ATLES  model  (Yoimg  1977)  and  unpublished  data 
(Manning  1985). 

The  transmittance  through  the  gaseous  effluents  of  the  fire  plume  and  the 
molecular  contribution  to  the  emitted  radiance  are  calculated  by  a  band  model 
adapted  from  the  ATLES  code  (Manning  1985).  For  single  wavelength 
calculations,  the  waveband  is  extremely  narrow.  The  spectral  resolution  is  not 
high  but  is  dependent  on  the  band  model  data  base  used  and,  generally,  is 
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25  cm’'.  An  important  feature  of  the  ATLES  band  model  is  its  ability  to  treat 
temperature  dependence  if  the  data  base  contains  parameters  for  high 
temperatures.  This  is  the  primary  reason  for  choosing  the  ATLES  model  rather 
than  the  LOWTRAN  code  (Kneizys  et  al.  1983). 

The  band  model  was  extracted  from  the  ATLES  code  and  adapted  for  use  in 
FITTE.  Code  to  perform  calculations  through  aircraft  and  missile  plumes  was 
removed,  large  arrays  were  reduced  in  size  by  limiting  the  number  of  species 
the  code  can  handle  simultaneously  and  by  limiting  the  size  of  the  data  base 
files  that  can  be  loaded,  and  the  parts  of  the  code  that  handle  Doppler  and 
Voigt  line  shapes  were  eliminated  because  only  Lorentz  line  shapes  are  needed 
for  FITTE. 

The  importance  of  temperature  dependence  was  described  (Manning  1985). 
Manning  points  out  that,  at  high  temperatures,  transmittance  is  not  strongly 
affected  in  the  2500-cm''  region,  but  it  is  greatly  reduced  in  the  2000-  and 
3300-cm’'  regions  as  hot-band  transitions  from  nearby  carbon  dioxide  and  water 
vapor  bands  become  active.  Because  the  radiance  emitted  in  a  spectral  band 
increases  as  transmittance  decreases,  correct  prediction  of  transmittance  in  the 
fire  plume  is  important. 

The  band  model  parameters  included  for  use  in  FITTE  are  taken  from  four  files: 
two  for  carbon  dioxide  and  two  for  water  vapor.  The  files  were  combined  into 
a  single  file  called  BPARAM.DAT.  The  carbon  dioxide  section  of  the  file 
covers  the  spectral  region  2000  to  2400  cm  ',  which  includes  the  important 
carbon  dioxide  fundamental  band  at  4.3  pm  (2300  cm’')  and  the  region  from 
3100  to  3700  cm  '.  It  includes  temperature  data  at  10  temperatures  from  100 
to  3000  K.  The  water  vapor  section  of  the  file  covers  the  spectral  regions  from 
50  to  2500  cm’'  and  from  2500  to  4500  cm’'.  It  contains  data  at  seven 
temperatures  between  300  and  3000  K.  The  BPARAM.DAT  file  represents 
most  of  the  important  absorption  bands  of  carbon  dioxide  and  water  vapor  in 
the  IR.  The  carbon  dioxide  files  were  supplied  with  the  ATLES  code.  The 
water  vapor  files  are  later,  unpublished  data  (Manning  1985). 

When  the  input  file  specifies  that  calculations  are  to  be  performed  for  a  spectral 
band,  propagation  is  calculated  using  ATLESF  for  a  number  of  wave  numbers 
across  the  band.  The  values  are  used  to  predict  values  for  a  selected  sensor 


response  function  for  the  waveband.  The  radiance  is  integrated  over  the  band, 
and  the  transmittance  is  averaged  over  the  band.  Calculations  are  made  using 
equations  (27)  and  (28). 

If  I(v)  represents  spectral  radiance  (W/mVsr/cm),  and  f(v)  represents  an 
instrument  spectral  response  function  with  a  range  of  zero  to  one,  the  integrated 
radiance  is  given  by 


/  =  /  /(V)  /(V)  dv 


(27) 


and  has  units  of  W/mVsr. 

The  average  transmittance  over  the  band  is 


f  ®(v)  /(v)  dv 

0  =z  1 - 

/  /(v)  dv 


(28) 


Note  the  difference  in  the  two  calculations.  The  radiance  is  an  integrated 
quantity;  if  the  radiance  is  constant  with  wave  number  and  the  instrument 
response  function  is  a  rectangular  function,  increasing  the  size  of  the  rectangular 
window  should  increase  the  radiance.  The  transmittance  is  an  average  over  the 
band;  if  the  transmittance  is  constant  with  wave  number  and  the  size  of  the 
rectangular  window  is  increased,  the  average  transmittance  should  not  change. 

2.2.3  Fire  Plume  Turbulence  Effects  on  Laser  Beam  Propagation 

Prediction  of  the  effects  of  fire  plume  turbulence  on  propagation  requires  the 
use  of  a  stochastic  model.  The  turbulent  hot  plume  can  be  regarded  as  a 
variable  density  lens  that  changes  with  time.  The  driving  forces  for  the 
fluctuations  are  entrainment  and  turbulent  mixing  of  ambient  temperature  air 
into  the  plume.  It  is  assumed  that  existing  theory  for  laser  beam  propagation 
through  the  turbulent  atmosphere  can  be  used  if  appropriate  changes  are  made 
to  the  definitions  of  parameters  and  the  ranges  in  which  they  are  valid.  The 
key  assumption  is  that  the  index  of  refraction  structure  constant  Cl  must  be 
replaced  by  a  structure  function  within  the  plume  dependent  on  the 
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time-averaged  temperature  distribution.  The  gradient  of  the  time-averaged 
temperature  distribution  determines  the  step  size  that  must  be  used  to  obtain 
accurate  results  from  the  calculations. 

The  turbulence  effects  model  predicts  values  of  laser  beam  parameters  at  the 
target  that  depend  on  the  turbulence-induced  fluctuations  in  the  LOS  index  of 
refraction.  It  is  a  modification  of  a  model  developed  (Thompson  and 
DeVore  1981)  for  ASL  based  on  work  in  the  literature  (Fante  1975).  The 
magnitudes  of  the  propagation  effects  are  a  complex  function  of  the  wavelength 
and  of  the  relative  sizes  of  the  laser  beam  and  the  plume. 

The  index  of  refraction  structure  constant  (Tatarski  1971)  is  given  by 


(29) 


where  n(r)  is  the  instantaneous  fluctuation  in  index  of  refraction  at  point  r-^  and 
ri2  is  the  distance  between  points  r,  and  ^2- 


The  brackets  denote  an  ensemble  average  over  all  points  (in  a  localized  region) 
separated  by  the  distance  r^2  •  In  practice,  an  ergodic  hypothesis  is  invoked  to 
replace  the  ensemble  average  by  a  time  average. 

The  index  of  refraction  of  air  is  given  by  (Strohbehn  1978) 


n  =  1.0  +  77.6x10'^ 


P 

T 


,  .  7.52x10-3 

1.0  +  - 


(30) 


where  P  is  the  pressure  (mbar),  T  is  temperature  (K),  and  A.  is  the  wavelength 
of  interest  (pm). 

A  humidity-dependent  term  negligible  at  visual  and  IR  wavelengths  was 
omitted.  Examination  of  equation  (30)  shows  that  temperature  fluctuations 
cause  index  of  refraction  fluctuations. 

For  atmospheric  turbulence  in  which  the  temperature  fluctuations  are  of  the 
order  of  1  K,  the  variance  of  the  index  of  refraction  is  proportional  to  the 


variance  of  the  temperature.  For  turbulence  in  fire  plumes  in  which 
temperature  fluctuations  are  of  the  order  of  tens  or  himdreds  of  degrees,  the 
variance  of  the  index  of  refraction  must  be  related  to  that  of  the  reciprocal  of 
the  temperature; 

where 

-  l)T. 


(31) 

(32) 


The  index  of  refraction  structure  constant  is  obtained  from  the  relation 


C/  =  1.91 


(33) 


where  is  the  outer  scale  of  turbulence  (Tatarski  1971). 


Effects  of  turbulence  on  propagation  in  the  atmosphere  are  calculated  using  the 
assumption  of  Kolmogorov  spectral  distribution  of  index  of  refraction 
fluctuations.  is  the  outer  scale  of  turbulence  (the  size  of  the  largest  eddies 
for  which  the  energy  cascade  to  smaller  eddies  is  linear),  and  is  the  inner 
scale.  For  the  atmosphere,  is  of  the  order  of  tens  of  meters  and  /„  is  of  the 
order  of  millimeters.  is  actually  an  inverse  function  of  height  above  the 
ground  for  neutral  and  unstable  Pasquill  stability  categories  and  is 
approximately  constant  for  a  stable  atmosphere  (Wyngaard  1973). 

For  turbulence  in  fire  plumes,  the  Kolmogorov  spectral  distribution  is  used  with 
the  outer  scale  set  equal  to  the  radius  of  the  plume.  Values  of  Cl  inside  the 
warm  region  of  the  plume  are  typically  one  to  three  orders  of  magnitude  higher 
than  ambient  values.  To  evaluate  Cl  along  the  path,  it  is  necessary  to  assume 
a  temperature  probability  distribution  about  the  mean  temperature  at  any  given 
point.  The  time-averaged  temperature  distribution  in  the  plume,  in  a  plane 
perpendicular  to  the  centerline,  is  a  circularly  symmetrical  Gaussian  distribution 
with  the  following  form: 


T  =  (Tc  -  T^)  ^ 


(34) 
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where  T  is  the  time-averaged  temperature  at  point  p,  Tc  is  the  centerline 
temperature,  is  the  ambient  temperature,  r  is  the  distance  of  point  p  from  the 
centerline,  and  <3(D)  is  the  diffusion  function  that  specifies  the  local  plume 
width. 

The  centerline  temperature  is  obtained  from  the  plume  model,  and  the  average 
temperature  at  point  p  is  used  to  predict  the  value  of  Cl(p)  in  the  turbulence 
effects  calculation.  In  actual  calculations,  the  simplified  exponential  expression 
of  equation  (33)  is  replaced  by  the  exponential  terms  of  equation  (17). 

FITTE  uses  established  turbulence  theory  (Fante  1975)  to  find  the  diffraction 
limited  spot  radius  (p^),  lateral  coherence  length  (pj,  short-  (pj  and  long-time 
(Pl)  average  spot  radii,  and  average  instantaneous  distance  of  the  centroid  (p^.) 
from  its  long-time  average  position.  The  parameters  are  used  to  predict  a 
centroid  jitter  radius  and  a  spot  area  magnification  v|/. 

The  increase  in  the  instantaneous  beam  area  at  the  target  is  calculated  from 

(  P 

111  -  ^  (35) 

where  is  the  diffraction  limited  spot  radius  at  the  target. 

This  result  can  be  used  to  estimate  the  average  decrease  in  intensity  at  the 
target.  No  detailed  information  about  spot  intensity  distributions  is  predicted 
because  the  distributions  are  a  function  of  relative  beam  and  plume  sizes,  the 
variations  in  path  integrated  particulate  concentration  for  incremental  areas  of 
the  spot  and  the  fluctuating  lens  represented  by  the  plume. 

For  the  points  outside  the  plume,  the  simplifying  assumption  made  is  that 
is  constant  and  equal  to  1.0  x  lO''"*  m'^^.  The  model  is  not  sensitive  to  the 
value  chosen. 

For  plume  temperatures  less  than  10  °C  above  ambient,  an  approximation  is 
used  to  evaluate  the  variance  of  the  index  of  refraction.  For  equation  (36) 
(George,  Alpert,  and  Tamanini  1977) 


a/  =  0.4  (f  -  T,) 


(36) 


where  T  is  the  mean  temperature  at  a  point, 
and  the  approximate  relation 


a 


2 

n 


(37) 


can  be  used  in  equation  (33)  to  find  C^.  The  values  obtained  for  al  using  this 
approximation  give  satisfactory  results  in  the  prediction  of  centroid  jitter  when 
compared  with  experiments  (D.  Bruce,  unpublished  data). 

When  the  turbulence  is  not  yet  folly  developed  within  a  distance  from  the  fire 
of  four  times  the  source  radius,  a  separate  method  is  used  to  calculate  C^.  The 
method,  described  elsewhere  (Thompson  and  Devore  1981),  relates  the  variance 
of  the  temperature  structure  function  to  the  mean  temperature  gradients  in 
the  plume  and  assumes  partial  mixing  of  ambient  air  in  this  region.  Cf  is 
calculated  from  equation  (33),  but  equations  (36)  and  (37)  are  replaced  by 


a/(A  r)  =  31, 


2  (  vr(A  r)  y 
’  1.91 


(38) 


and 


a 


2 

n 


(39) 


where  T/^  is  the  temperature  of  the  partially  mixed  fractional  volume  of  air. 

This  method  of  calculation  is  invoked  automatically  by  the  model  when  the 
LOS  passes  through  the  incompletely  mixed  region  of  the  plume. 
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2.3  Model  Implementation 
2.3.1  FITTE 

FITTE  performs  calculations  to  predict  effects  for  specific  LOS. 
FITTE/FGLOW  has  three  defined  scenarios.  Scenarios  1  and  2  are  single  LOS 
(FITTE)  scenarios  used  to  predict  the  effects  of  the  fire  plume  on  general 
systems  such  as  laser  designators,  and  laser-  or  other  seeker  systems. 
Scenario  3,  discussed  below  as  FGLOW,  performs  calculations  for  multiple 
LOSs  to  predict  radiant  images  of  fires  that  may  affect  imaging  systems  such 
as  IR  cameras. 

Scenarios  1  and  2  have  identical  geometry,  as  seen  in  figure  5,  but  differ  in  the 
direction  of  the  turbulence  calculation.  Separate  calculations  are  needed 
because  calculation  results  depend  on  path-weighted  averages.  The  target  and 
emitted  radiance  are  calculated  as  seen  from  the  observer  position  for  both 
scenarios.  For  scenario  1,  the  turbulence  effects  are  calculated  at  the  target 
position;  for  a  laser  beam  propagating  from  observer  to  target.  For  scenario  2, 
they  are  calculated  at  the  observer  position;  for  a  reflected  laser  beam 
propagating  from  target  to  observer. 

2. 3. 1.1  Scenario  1. — Scenario  1  corresponds  most  nearly  to  the  original  version  of 

FITTE,  and  it  is  the  default  for  running  the  model.  The  observer  can  be  any 
of  the  following:  a  laser  designator,  an  observer  with  night-vision  goggles,  or 
a  visible  or  IR  imaging  system.  The  target  has  no  assigned  physical  dimensions 
but  may  have  a  nonambient  temperature  and  emissivity.  The  radiative  transfer 
calculations  are  performed  so  the  values  produced  are  given  at  the  observer 
position.  The  transmittance  applies  at  either  end  of  the  path.  Unless  a 
waveband  calculation  is  specified,  turbulence  effects  parameters  for  the  beam 
(assumed  to  be  a  laser  beam)  are  predicted  at  the  target  position. 

Various  combinations  of  output  parameters  can  be  used  to  predict  quantities 
such  as  the  relative  laser  intensity  on  the  target  (initial  intensity  times 
transmittance  divided  by  spot  magnification),  target  contrast  (target  radiance 
divided  by  the  sum  of  target  and  emitted  radiances),  etc. 
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Figure  5.  Geometry  for  the  single  LOS  scenarios.  The  helicopter  carries  a  laser 
designator  for  scenario  1  or  a  laser-guided  munition  for  scenario  2. 

2.3. 1.2  Scenario  2. — Scenario  2  is  used  much  the  same  way  as  scenario  1  and  performs 

the  same  calculations  for  radiative  transfer.  It  was  included  for  the  calculation 
of  turbulence  effects  on  a  laser  beam  reflected  fi-om  the  target  to  the  observer. 
The  effects  are  different  from  the  effects  in  scenario  1  for  all  cases  in  which  the 
fire  plume  is  not  at  the  center  of  the  propagation  path,  because  the  turbulence 
in  the  plume  distorts  the  wave  front  in  a  manner  similar  to  a  set  of  diverging 
lenses;  therefore,  the  net  effects  depend  on  the  path-weighted  plume  turbulence. 

Time-series  calculations  can  be  performed  for  these  scenarios  to  give  tabular 
output  of  the  parameters.  Figure  6  shows  graphs  of  the  predicted  fluctuations 
of  several  parameters. 
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Figure  6.  Graphs  of  model  output  data  for  a  time-series  calculation  at  10.6  fim.  The 
dashed  lines  show  the  values  predicted  by  the  mean-value  model. 


36 


2. 3. 1.3  Geometric  Relations  for  Plume  Parameters  on  the  LOS. — ^Because  integration 

along  the  LOS  is  required  to  evaluate  the  fire  or  plume  effects,  geometrical 
relations  are  needed  to  predict  plume  parameter  values  at  points  on  the  LOS. 
The  plume  parameters  are  expressed  as  functions  of  the  distance  along  the 
centerline  and  of  the  radial  distance  from  the  centerline  to  the  point  on  the 
LOS. 

Any  point  p  on  the  LOS  has  a  projection  (x-^j  0,  zj  on  the  plane  of  the  time- 
averaged  plume  centerline,  and  the  shortest  distance  from  the  centerline  to  the 
projected  point  is  along  a  line  perpendicular  to  the  centerline.  The  equation  of 
any  such  line  has  the  form 


z  =  mx  +  b 


(40) 


where  m  is  the  negative  inverse  of  the  centerline  slope  at  a  point  (x^,  0,  Zcf, 


m  = 


1 


dzidx  j 


(41) 


Because  (x^  0,  zf)  is  also  on  the  line,  the  intercept  h  is  given  by 


h  = 


3 

2 


(42) 


An  equation  for  Zq,  the  centerline  point  closest  to  (x-^  0,  z^,  is  found  by 
substituting  equations  (41)  and  (42)  into  equation  (40)  and  replacing  Xc,  by 
substitution  from  the  centerline  equation.  The  equation  obtained  is 


3 

2 


Zy  Zi 


=  0.  (43) 


Equation  (43)  can  be  solved  by  the  substitution 

Z  =  ^  Zqj  +  Zy 


(44) 
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and  subsequent  use  of  the  solution  for  quartic  equations  found  in  the  CRC® 
Standard  Mathematical  Tables  (Beyer  1981). 

The  distance  of  any  point  on  the  LOS  from  the  closest  centerline  point  can  be 
found  and  used  to  evaluate  the  plume  parameters. 

The  distance  from  the  centerline  point  to  the  projection  of  the  LOS  point  on  the 
x-z  (wind-aligned)  plane  is  denoted  by  h.  The  distance  from  the  plume 
centerline  to  the  ground  along  the  same  line  is  defined  as  h^.  When  is 
greater  than  2.0  a,  the  constants  and  Kj  in  the  plume  equations  are  set  to  the 
values  for  a  vertical  plume  because  no  reflection  from  the  ground  occurs. 

2. 3. 1.4  Information  about  Calculations  along  the  LOS. — ^For  radiance  and  transmittance 

calculations,  the  default  number  of  calculation  intervals  through  the  plume  is 
30,  giving  a  step  length  of  0.2  g.  This  does  not  affect  the  model  accuracy  for 
most  LOS.  Step  size  effects  are  most  likely  to  occur  for  LOS  near  the  plume 
centerline  and  very  close  to  the  fire.  The  calculation  of  turbulence  effects  is 
performed  using  an  adaptive  Simpson’s  rule  integration  subroutine,  ASRI, 
rather  than  a  fixed  number  of  calculation  intervals.  A  single  calculation  interval 
is  used  between  the  plume  and  the  observer  when  molecular  absorption 
calculations  are  performed. 

When  waveband  calculations  are  performed,  the  molecular  effects  subroutines 
(the  ATLESF  set)  are  called  for  the  first  interval,  and  subsequent  calls  are  made 
only  when  the  temperature  change  from  the  previous  interval  exceeds  a  fraction 
(default  =  0.1)  of  the  previous  value. 

2,3.2  FGLOW:  Fire/Plume  Radiant  Emission  and  Transmittance 

FGLOW  is  an  extension  of  FITTE  that  performs  calculations  for  multiple  LOSs 
to  generate  a  simulated  image  matrix.  It  does  not  calculate  turbulence  effects 
on  laser  propagation. 

Figure  7  shows  scenario  3,  the  FGLOW  scenario.  An  imaging  system  is 
located  at  the  observer  position  with  its  LOS  specified  by  a  target  position,  and 
a  fire  or  a  plume  segment  is  within  the  field  of  view.  FITTE  calculates  the 
emitted  radiance  for  an  array  of  pseudotarget  positions  to  create  a  file  that 
contains  a  matrix  representation  of  the  radiant  image  of  the  fire  and 
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transmittance  of  the  plume.  The  array  is  defined  by  an  input  card  to  have 
n  rows,  m  columns,  and  an  angular  resolution  (in  milliradians)  between 
elements. 


Figure  7.  Top  view  of  scenario  3 
(FGLOW)  geometry.  Situation  shown  is 
for  the  240'’  wind  direction  of  example  9. 


Because  the  calculation  matrix  for  FGLOW  typically  contains  many  LOSs  that 
do  not  pass  throu^  the  plume,  the  model  sets  a  flag  and  stores  the  values 
calculated  for  the  first  LOS  that  does  not  intercut  the  plume.  Calculation 
speed  is  increased  by  using  these  stored  values  for  all  subsequent  LOS  that  do 
not  intercept  the  plume. 
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3.  Caveats 


3.1  Grade  of  Software 

The  FITTE  model  should  be  classified  as  a  research  model.  This  choice  of 
software  grade  is  based  on  the  fact  that  no  quantitative  evaluations  were  made 
for  the  spatial  and  time  variations  introduced  into  the  model  and  for  the  limited 
evaluations  of  predictions  of  beam  wander  and  spread. 

3.2  Model  Failure 

The  model  fails  to  run  if  selected  input  values  do  not  meet  acceptance  criteria. 
Failure  occurs  if:  (1)  the  observer  and  target  positions  are  on  the  same  side 
of  the  wind-aligned  axis  in  the  model  coordinate  system,  or  (2)  either 
waveband  limit  on  the  BAND  card  is  outside  the  range  of  the  molecular 
absorption  data  in  BPARAM.DAT. 

The  model  fails  to  run  if;  (1)  the  files  BPARAM.DAT  and  FIT4D.DAT  are 
not  located  in  the  same  directory  as  the  executable  program,  (2)  a  non- 
EOSAEL  format  input  card  is  encountered,  or  (3)  an  unexpected  end-of-file  is 
encountered  in  an  input  file. 

The  model  fails  if  the  disk  runs  out  of  space  for  the  output  files.  The  model 
is  most  likely  to  fail  if  FGLOW  is  run  for  a  large  image  matrix  or  a  time-series 
of  smaller  FGLOW  images. 

Caution  must  be  exercised  under  the  following  conditions: 

•  The  model  runs  if  many  of  the  possible  input  parameters  are  out  of  bounds. 
In  this  case,  a  default  value  is  substituted  for  the  out-of-bound  parameter,  the 
output  file  contains  a  message  stating  DEFAULT(S)  REPLACE  INVALID 
INPUT  ON  CARD  CARDNAME,  and  the  default  values  appear  in  the  printout 
of  input  parameters  contained  in  the  output  file. 

•  Bounding  values  are  set  for  a  parameter  if  the  value  is  known  to  be 
physically  bounded  (target  emissivity  must  be  between  0  and  1)  or  if  an 
out-of-bound  value  affects  the  accuracy  of  model  computations  (windspeed  must 
be  greater  than  0.10  because  it  appears  in  the  denominator  of  several  equations). 
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•  The  model  performance  was  tested  only  under  relatively  benign  atmospheric 
conditions:  windspeeds  less  than  about  10  m/s  (22  mph)  on  clear  or  cloudy  days 
with  unstable  or  neutral  stability.  Model  performance  is  expected  to  degrade 
under  more  severe  weather  conditions,  but  it  is  not  known  how  rapid  or  severe 
the  degradation  will  be.  No  provision  exists  for  including  effects  of  obscurants 
such  as  rain,  fog,  or  smoke  clouds  with  the  effects  of  the  fire  plume. 

3.3  Verification  Tests 

Model  verification  was  obtained  from  laboratory  and  field  tests.  All 
comparisons  to  date  were  made  with  the  predictions  of  the  Gaussian  plume 
model.  This  is  the  most  meaningful  approach  to  the  question  of  general  model 
validity  because  time-averaged  data  correspond  well  with  a  Gaussian  pluirte 
(D.  Bruce,  unpublished  measurements).  Transmittance,  radiance,  and  plume 
centerline  shape  and  spread  for  the  time-averaged  plume  agree  well  with  field 
test  data.  The  laser  wander  and  beam  spread  predictions  were  verified  by 
laboratory  experiments.  Aerosol-specific  extinction,  generally,  agrees  well  with 
laboratory  and  field  data,  but  visible  and  near-IR  values  of  albedo  are  low 
compared  to  measured  values. 

Comparisons  with  field  data  suffer  from  various  problems.  If  point 
measurements  are  made,  the  results  depend  strongly  on  plume  meander  because 
of  changing  wind  conditions.  If  the  instruments  are  at  fixed  positions,  it  is 
possible  to  obtain  little  or  no  data.  If  data  are  obtained,  auxiliary  measurements 
are  needed  to  determine  the  location  of  the  probe  within  the  plume  relative  to 
the  instruments.  If,  as  at  the  BICT  III  test,  a  crane  is  used  to  position  an 
instrument  cage  within  the  plume,  data  are  obtained,  but  the  data  reduction 
problem  of  relative  position  is  even  more  difficult.  Similar  problems  occur  for 
transmittance  measurements  that  represent  LOS  integrated  data.  As  a  result, 
limited  data  are  available  because  of  difficulties  encountered  in  making 
measurements  through  plumes  in  the  field  and  in  analyses  of  the  data. 

3.3.1  Fire  and  Plume  Properties 

The  shape  of  the  time-averaged  plume  centerline  (Briggs  1969)  was  based  on 
extensive  observations  of  plumes  from  smokestacks  under  a  variety  of 
meteorological  conditions.  The  applicability  of  this  model  to  ground-level 
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small  fires  was  verified  by  comparing  visible  plume  images  (photographs  or 
video)  with  model  predictions  for  appropriate  fire  size  and  wind  conditions. 
Good  agreement  was  reported  (Bruce  and  Sutherland  1986). 

The  plume  centerline  angle  of  elevation  depends  on  the  ratio  of  the  vertical 
velocity  to  the  (horizontal)  windspeed.  Plume  angles  of  elevation  as  a  function 
of  windspeed  were  reported  (Hall  and  Manning  1986)  from  imagery  taken  at 
BICT  III.  Digitized  images  from  each  of  10  data  segments  were  used  in  the 
analyses.  The  data  segments  were  chosen  to  represent  a  range  of  windspeeds. 
After  reclassifying  the  windspeed  for  two  segments  of  the  data,  the  authors 
found  a  linear  relation  between  angle  of  elevation  and  windspeed  for 
windspeeds  ranging  from  1  to  7  m/s.  The  authors  mention  that  uncertainties  are 
introduced  in  determining  the  position  of  the  origin  of  the  plume  and  the  need 
to  ignore  centerline  curvature  (Hall  and  Manning  1986). 

Figure  8  shows  the  measured  values  compared  with  model  predictions  for  the 
fire  type  (2)  closest  to  the  BICT  III  fires.  The  error  bars  indicate  the  standard 
deviation  of  the  data.  Agreement  is  reasonable,  but  the  difference  in  slope 
indicates  that  changes  in  the  model  would  be  appropriate.  The  changes  should 
probably  involve  an  increase  in  thermal  flux  as  a  function  of  windspeed. 
Experiments  are  currently  being  conducted  to  determine  effects  of  windspeed 
on  fire  and  plume  parameters. 

Another  important  measure  for  the  plume  model  is  how  well  it  predicts  plume 
spread.  The  ability  of  the  model  to  predict  plume  spread  can  be  verified  by 
comparisons  of  total  angular  spread  observed  in  images  with  the  value  obtained 
from  the  entrainment  parameter  ^  for  various  windspeeds.  This  type  of 
comparison  was  made  for  a  number  of  plume  images,  and,  with  the  assumption 
that  the  visible  plume  edge  is  at  the  2c  point  of  the  model  Gaussian 
distribution,  agreement  is  good.  This  means  that  the  model  predictions  are 
typically  within  the  range  of  uncertainty  for  the  measured  values. 
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Figure  8.  Comparison  of  measured  and 
modeled  plume  elevation  angles. 


3.3.2  Aerosol  Optical  Properties 

Laboratory  and  field  measurements  of  aerosol  optical  properties  in  plumes  from 
diesel  fiiel  fires  were  reported.  Figure  9  shows  a  comparison  of  extinction 
coefficients  from  FITTE  with  laboratory  data  (Bruce  et  al.  1991).  The 
measurements  were  made  with  a  flowthrough  photoacoustical  system  with  an 
incorporated  short-path  transmissometer  from  which  data  is  obtained  for  both 
specific  absorption  and  extinction  coefficients.  Correlated  measurements  of 
aerosol  characteristics  were  obtained  by  timed  collection  of  filter  samples  and 
subsequent  analyses  to  obtain  total  particle  weight,  particle  size  distribution,  and 
chemical  identification.  Error  bars  indicate  variations  of  measured  values  with 
the  naturally  occurring  variation  of  aerosol  size  distribution  during  a  series  of 
experiments.  Agreement  at  10.6  pm  is  due  to  the  fact  that  the  measured  values 
at  that  wavelength  were  used  to  determine  the  absolute  scale  for  the  Mie 
calculations  used  for  the  model  optical  parameters  at  other  wavelengths. 
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modeled  extinction  coefficients. 

Two  series  of  field  measurements  (Bruce  and  Richardson  1983; 
Bruce  et  al.  1989)  of  optical  properties  of  diesel  fire  plume  aerosols  at  10.6  pm 
gave  extinction  coefficients  averaging  about  20  percent  higher  than  those  for  the 
laboratory  measurements  and  showing  a  greater  variability  than  the  laboratory 
data.  We  attribute  the  increased  values  to  the  presence  of  a  greater  number  of 
large  aerosols  in  the  measured  size  distributions.  The  values  from  the  field  data 
(1.3  +  0.7  and  1.2  +  0.3  mVg)  overlap  the  laboratory  value  (0.99  ±  0.12  m^/g). 

The  modeled  and  measured  values  of  single  scattering  albedo  are  compared  in 
table  3.  The  measured  values  are  from  the  same  laboratory  experiments  as  the 
specific  extinction  and  are  calculated  using  the  formula  co  =  1  -  a/s.  No  value 
was  obtained  at  1.06  pm  because  no  gas  that  absorbs  strongly  at  this 
wavelength  is  available  to  provide  calibration  data  for  the  spectrophone,  and 
changes  in  beam  geometry  for  different  wavelengths  preclude  the  use  of 
calibration  from  other  wavelengths  (C.  W.  Bruce,  private  communication). 
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Table  3.  Comparison  of  albedo  values 


X  (|im) 

’^MODEL 

^MEASURED 

0.488 

0.25062 

0.34720 

3.390 

0.22174 

0.32540 

10.600 

0.15000 

0.15000 

The  substantial  difference  between  the  modeled  and  measured  values  at  shorter 
wavelengths  indicates  problems  with  the  values  currently  used  in  the  model. 
Revision  of  the  values  await  data  at  additional  wavelengths,  since  the  functional 
dependence  is  highly  non-linear  in  the  short  wavelength  region.  The  use  of  Mie 
calculations,  valid  only  for  spherical  particles,  and  particle-to-void  theory  for 
the  stringy  aerosol  particles  is  subject  to  large  uncertainties,  as  was  the 
assumption  that  the  particles  could  be  described  by  a  bimodal  size  distribution. 
The  primary  effect  on  model  output  is  a  change  in  predicted  emitted  radiance 
at  visual  and  near-IR  wavelengths.  Experimental  runs  using  measured  values, 
and  with  scaled  increases  in  albedo  values  at  intermediate  wavelengths,  indicate 
that  emitted  radiance  decreases  by  12  to  15  percent  in  the  visual  and  near-IR 
wavelength  regions. 

3.3.3  Plume  Distribution  Functions 

Radiant  emission,  measured  with  thermal  imaging  systems  during  the  Battlefield 
Emissive  Sources  Trial,  Oldebruck,  the  Netherlands,  Europe  (BEST  ONE)  test, 
was  used  to  check  the  predictions  of  FGLOW.  Images  were  analyzed  to  give 
the  time  dependence  of  the  total  area  that  exceeded  various  blackbody 
equivalent  temperatures  (S.  B.  Crow,  private  communication).  Good  agreement 
of  measured  and  simulated  radiation  temperature  spatial  distributions  was 
obtained  for  top  hat  source  temperatures  ranging  from  1223  to  1473  K  (Bruce 
1989).  This  is  roughly  300  K  higher  than  the  value  currently  used  for  the 
model.  Peak  temperatures  measured  close  to  the  edge  of  one  of  the  fuel  drums 
during  the  BICT  III  test  were  approximately  1250  K  (D.  Bruce,  unpublished 
data).  Calibrated  imagery  taken  at  BEST  ONE  with  short-  and  long-wavelength 
cameras  and  several  narrow-band  filters  (Leidner  and  Clement  1987)  show 
maximum  radiation  temperatures  between  800  and  1300  °C  for  an  open  tray  fire 
and  between  600  and  900  °C  for  burning  vehicles  of  which  the  flames  were 
mostly  enclosed. 
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Peak  vertical  plume  velocities  measured  at  the  BICT  III  test  (Bruce  et  al.  1989) 
were  approximately  2.9  m/s  at  3  m  from  the  fire  and  1.2  m/s  at  approximately 
14  m  from  the  fire. 

Molecular  combustion  product  concentrations  in  diesel  fire  plumes  were 
measured  (Bruce  et  al.  1989;  Kaaijk  1986).  Maximum  measured  carbon 
dioxide  concentrations  ranged  from  100  ppm  to  >1000  ppm  above  the  ambient 
concentration,  and  the  maximum  measured  water  vapor  increase  ranged  from 
0  to  7  mb  above  ambient.  Other  effluent  gases  occur  in  very  low 
concentrations.  Benzene  was  the  prevalent  aromatic  hydrocarbon,  with  a 
measured  concentration  range  of  0.3  to  3.0  mg/m^  (Kaaijk  1986),  and  carbon 
monoxide  concentration  was  less  than  100  ppm  because  it  was  not  detected  by 
the  chemical  analyses.  The  concentrations  reported  by  Kaaijk  represent  an 
average  value  for  the  plume  position  during  the  sampling  time,  while  those 
reported  by  Bruce  et  al.  include  instantaneous  as  well  as  time-averaged  values. 
The  concentrations  in  both  cases  were  measured  approximately  15  m  from  the 
fire.  The  measurements  can  only  be  used  to  determine  if  the  effluent 
concentrations  used  in  FITTE  are  approximately  correct. 

The  carbon  dioxide  concentration  was  chosen  for  comparison  with  the  model 
because  it  is  a  more  direct  measurement.  The  assumption  was  made  that  the 
mean  concentration  at  a  distance  of  15  m  from  the  base  of  the  plume  was 
approximately  667  ppm.  Data  given  for  plume  spread  and  fire  size 
(Bruce  et  al.  1989)  were  used  to  estimate  values  at  the  source  and  obtain  a 
mean  concentration  of  approximately  8000  ppm;  much  less  than  the 
136,000  ppm  modeled  in  earlier  versions  of  FITTE.  The  very  large  number 
was  obtained  based  on  the  assumption  (Manning  1985)  that  the  effluent  gases 
represented  100  percent  of  the  gases  above  the  fire  and  the  fact  that  the  effluent 
gas  concentrations  were  not  modeled  as  dependent  on  the  bum  rate.  The 
current  model  for  effluent  gas  concentrations  ties  their  values  to  the  rate  of  fuel 
consumption  and  results  in  source  concentrations  of  carbon  dioxide  between 
6000  and  17,000  ppm,  in  satisfactory  agreement  with  the  measured  values. 

The  aerosol  mass  concentrations  reported  from  filter  data  (Kaaijk  1986)  ranged 
from  0.001  to  0.036  g/m^  and  were  not  corrected  for  the  time  the  sampler  was 
not  in  the  plume  or  for  the  fact  that  the  sampling  system  was  usually  near  the 
edge  of  the  plume  when  it  was  within  the  plume.  The  time  in  the  plume  varied 
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from  less  than  10  to  70  percent  of  the  20-min  sampling  period.  The  peak 
aerosol  mass  concentrations  derived  from  optical  measurements  within  the 
plume  (Bruce  et  al.  1989)  were  typically  0.500  g/m^  with  an  observed 
maximum  value  greater  than  0.640  g/ml  The  same  estimation  procedure  used 
for  the  source  gaseous  effluent  concentration  gives  a  source  aerosol  mean  mass 
concentration  of  approximately  2.6  g/m^;  approximately  20  percent  higher  than 
the  value  used  in  the  model. 

A  less  direct  indicator  of  the  adequacy  of  the  plume  distribution  functions  is 
given  by  comparison  of  measured  and  simulated  transmittance.  The  comparison 
depends  on  accurate  optical  parameters  as  well  as  the  distribution  functions  and 
must  be  performed  with  appropriate  matching  of  beam  parameters. 
Comparisons  of  measured  transmittance  at  BEST  ONE  with  FITTE  predictions 
were  reported  (Gillespie  1987;  Bruce  1989).  The  results  for  this  test  are 
compromised  by  the  fact  that  no  survey  data  are  available  to  provide  reliable 
position  data. 

Gillespie  used  a  modified  version  of  FITTE  (Manning  1985)  available  in  early 
1987  with  calculations  from  the  PC  version  of  LOWTRAN  to  account  for 
transmittance  effects  of  natural  aerosols.  This  version  of  the  model  did  not 
perform  area  averaging  for  large  beam  sizes.  Trials  were  divided  into  two 
classes  according  to  the  magnitude  of  transmittance  change  during  the  trial. 
Results  of  comparing  simulated  single  wavelength  visible  with  measured 
broadband  visible  data  were  poor,  typically  differing  by  more  than  20  percent. 
Results  at  1.06  pm  were  generally  better,  differing  within  10  percent  for  five 
trials  and  greater  than  20  percent  for  four.  In  the  3-  to  5-pm  band,  differences 
were  less  than  10  percent  for  five  trials,  less  than  20  percent  for  three  trials,  and 
greater  than  20  percent  for  two  trials.  In  the  8-  to  14-pm  band,  differences 
were  less  than  10  percent  for  two  trials,  less  than  20  percent  for  three  trials,  and 
greater  than  20  percent  for  four  trials.  Gillespie  categorizes  the  results  as  a  fair 
validation  of  FITTE. 

A  number  of  factors  probably  account  for  some  of  the  disagreement.  The 
measured  transmittance  at  1 .06  pm  is  lower  than  in  the  visible  band  for  most 
of  the  tests  used  in  the  comparison.  A  later  test  (Farmer  et  al.  1989)  suggests 
that  the  results  could  be  caused  by  vignetting  of  the  beam  at  one  of  the 
wavelengths,  operation  of  a  detector  in  a  nonlinear  response  range,  or 


misalignment  of  the  optical  system.  For  all  wavelengths,  there  were  a  number 
of  trials  for  which  no  comparison  was  made  because  the  model  indicated  that 
the  LOS  did  not  pass  through  the  plume.  This  indicates  problems  with  the 
available  input  data  and  may  have  affected  all  aspects  of  the  comparison. 

Bruce  used  the  revised  version  of  the  1987  FITTE  module  to  perform 
comparisons  for  one  of  the  BEST  ONE  trials.  She  varied  the  fire  parameters 
to  obtain  good  agreement  with  calibrated  IR  imagery,  and  used  the  area 
averaging  feature  of  FGLOW  (the  angular  resolution  was  chosen  to  match  the 
transmissometer  beam  size)  to  find  transmittance  at  three  IR  wavelengths  for 
a  variety  of  positions  in  the  plume.  The  area  near  the  estimated  position  of  the 
transmissometer  LOS  was  examined,  and  an  area  was  found  in  which  agreement 
of  simulation  and  measurements  was  less  than  25  percent  for  all  wavelengths. 

3.3.4  Effects  of  Thermal  Plume  Turbulence  on  Laser  Beam  Propagation 

Concurrent  laboratory  measurements  were  made  of  temperature  fluctuations  in 
plumes  and  of  laser  beam  spread  and  jitter  for  a  beam  traversing  the  plume 
(Bruce  and  Sutherland  1986).  The  experiments  verified  the  method  used  to 
calculate  the  localized  turbulence  structure  function  within  the  plume.  The 
measured  plume  centerline  temperature  and  radius  were  used  to  tailor  the  FITTE 
plume  for  simulations  that  showed  that  the  FITTE  predictions  of  laser  beam 
wander  are  good  (within  10  percent)  and  predictions  of  beam  spread  are  good 
to  satisfactory  (within  15  to  25  percent)  for  regions  outside  the  flames.  The 
errors  are  greatest  in  the  region  just  outside  the  flames.  The  experimental  data 
identified  a  method  of  calculation  to  be  used  for  the  area  close  to  the  fire  in 
which  the  plume  is  not  yet  well  mixed. 
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4.  Operations  Guide 

FITTE  can  be  run  using  the  standard  EOSAEL  driver.  It  is  also  furnished  with 
a  dedicated  driver  for  use  on  a  PC.  If  the  driver  is  used,  link  the  FITTE 
subroutines  with  the  program  FITDRV.  FOR,  found  in  the  FIT  ACC  directory. 
FITDRV.FOR  contains  several  subroutines  and  a  block  data  file  EOSLBD. 
Several  example  auxiliary  programs  are  furnished  that  can  be  used  to  create 
pseudogreyscale  printouts  of  FGLOW  output  or  to  modify  the  FGLOW  output 
files  for  use  by  programs  that  can  create  contour  plots.  These  examples  are 
described  in  appendix  B. 

The  discussion  and  examples  in  the  manual  apply  to  the  code  and  information 
supplied  by  the  author  to  the  custodians  of  the  EOSAEL  library.  Some  of  them 
may  not  apply  to  the  actual  distribution.  Check  the  distribution  for 
documentation  files  that  may  give  information  on  changes  in  the  code.  The 
example  files  furnished  with  the  distribution  were  produced  with  that  version 
of  the  code. 

4.1  Input 

The  input  cards  (actually  input  file  records)  use  the  standard  EOSAEL  format 
(A4,  6X,  7E10.4).  Run  repeated  scenarios  with  the  usual  EOSAEL  procedure. 
Six  to  nine  data  cards,  normally,  are  used  with  GO  and  DONE  cards  controlling 
program  execution.  Tables  4  through  17  contain  descriptions  of  the  FITTE 
input  cards.  The  accompanying  text  gives  additional  information  about  some 
of  the  cards,  including  definitions  of  new  cards  and  changes  to  some  previously 
defined  cards. 

Running  the  model  with  the  standard  EOSAEL  driver  requires  an  EOSAEL 
WAVL  card  and  a  card  to  call  the  FITTE  module.  Usually,  it  also  requires  the 
cards  REFD,  SCEN,  SRCL,  METD,  and  DETD  or  SCN3.  SCEN  and  METD 
may  be  omitted  if  the  data  were  furnished  to  EOSAEL  main,  in  which  case, 
flags,  used  as  FITTE  calling  parameters,  indicate  the  need  to  convert  the  data 
units.  Use  all  other  cards  to  invoke  various  model  options. 
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REFD  contains  information  about  the  orientation  of  the  user  coordinate  system 
and  several  model  control  parameters.  The  position  data  for  the  observer  and 
target,  given  on  the  SCEN  card,  and  for  the  fire,  given  on  the  SRCL  card,  can 
be  in  any  convenient  right-handed  coordinate  system.  Give  the  x-axis  heading 
of  that  (user)  coordinate  system.  The  heading  and  the  wind  direction  on  the 
METD  card  are  given  by  angles  measured  clockwise  from  north.  Use  the 
control  parameters  on  REFD  as  indicated. 


Table  4.  The  REFD  card 


Identifier  Variable 

REFD 

Default 

Description 

Reference  data 

XHEAD 

x-axis  heading  of  user  coordinate  system 
(degrees  clockwise  from  north) 

ITYPE 

1 

code  for  fire  type  (1  through  4) 

1  -  jeep 

2  -  truck 

3  -  tank 

4  -  small  fuel  depot 

ISCN 

1 

code  to  specify  FITTE  scenario  (1  through  3) 

1  -  propagation  from  observer  to  target 

2  -  propagation  from  target  to  observer 

3  -  predict  propagation  to  imager 

lAVG 

0 

code  to  choose  time-averaged  plume 
calculation:  0  or  positive  number  =>  use  4D 
plume  negative  number  =>  use  time-averaged 
plume 

NTURB 

0 

negative  value  cancels  turbulence  calculations 
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The  SRCL  card  contains  the  position  of  the  origin  of  the  FITTE  coordinate 
system  in  the  user  coordinate  system.  The  x-axis  of  the  FITTE  coordinate 
system  is  aligned  with  the  wind  vector. 


Table  5.  The  SRCL  card 


Identifier  Variable 

Default 

Description 

SRCL 

Fire  location  data 

US(1) 

x-coordinate  of  center  of  fire  (m) 

US(2) 

y-coordinate  of  center  of  fire  (m) 

US(3) 

z-coordinate  of  base  of  fire  (m) 

The  SCEN  card  gives  the  observer  and  target  positions  in  the  user  coordinate 
system.  The  positions  must  lie  on  opposite  sides  of  the  x-axis  in  the  FITTE 
coordinate  system. 

Table  6.  The  SCEN  card 

Identifier  Variable  Default 

Description 

SCEN 

LOS  data 

UO(l) 

observer  x  coordinate  (m) 

UO(2) 

observer  y  coordinate  (m) 

UO(3) 

observer  z  coordinate  (m) 

UT(1) 

target  x  coordinate  (m) 

UT(2) 

target  y  coordinate  (m) 

UT(3) 

target  z  coordinate  (m) 

Most  of  the  input  data  for  the  METD  card  are  self-explanatory.  The  wind 
direction  is  determined  by  the  usual  meteorological  convention;  the  upwind 
direction  measured  clockwise  from  north.  The  Pasquill  stability  category  has 
a  range  of  1  to  6,  compared  with  the  usual  classification  1  =  A,  2  =  B,  etc. 
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Table  7.  The  METD  card 


Identifier  Variable 

METD 

Default 

Description 

Meteorological  data 

UBAR 

mean  ambient  windspeed  (m/s) 

WDIR 

wind  direction  (meteorological 
convention,  ®) 

TAIR 

ambient  air  temperature  (°C) 

RH 

relative  humidity  (%) 

RHO 

ambient  air  density  (g/m^) 

IPAS 

Pasquill  stability  category  (range 

1  to  6) 

The  DETD  card  specifies  the  initial  laser  heam  diameter  for  a  scenario  1  or  2 
simulation.  Use  any  convenient  diameter  because  the  model  predicts  beam 
magnification  rather  than  diameter  after  propagation. 


Table  8.  The  DETD  card 


Identifier 

Variable 

Default 

Description 

DETD 

Laser  data 

BEAM 

beam  diameter  (cm) 

Use  the  TARG  card  to  obtain  contrast  data  for  a  target  with  a  temperature  other 
than  ambient  for  a  scenario  1  or  2  situation. 


Table  9.  The  TARG  card 


Identifier  Variable 

TARG 

Default  Description 

Target  data  for  thermal  emission  calculations 

TTARG 

target  temperature  (°C) 

ETARG 

target  emissivity  (dimensionless, 

range  0  to  1) 
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Use  the  MOLS  card  to  specify  inclusion  of  molecular  effects  in  the  radiance 
and  transmittance  calculations.  It  also  controls  whether  the  transmittance  and 
radiance  calculations  are  performed  only  through  the  plume,  or  if  the  effect  of 
molecular  absorption  and  emission  in  the  atmosphere  between  plume  and 
observer  is  included.  The  ambient  carbon  dioxide  and  water  vapor 
concentrations  are  used  in  this  calculation;  other  atmospheric  gases  are  not 
included.  If  lATM  is  1,  the  stepwise  calculation  through  the  plume  is 
continued  from  the  near  plume  edge  to  the  observer  position. 


Table  10. 

The  MOLS  card 

Identifier 

MOLS 

Variable 

Default 

Description 

Molecular  effects  calculation  control 

IMOL 

1 

flag  for  molecular  transmittance  calculation: 
0  =>  molecular  transmittance  not  included 

1  =»  molecular  transmittance  included 

lATM 

0 

flag  for  calculations  outside  the  plume: 

0  =»  omit  calculations  between  observer  and 
plume 

1  =>  include  calculations  between  observer 
and  plume 

Use  the  BAND  card  to  specify  waveband  calculations  and  to  furnish  parameters 
for  their  control.  It  also  controls  the  instrument  response  function  used  for  a 
waveband  calculation.  There  are  three  instrument  response  ftinctions  from 
which  to  choose:  a  rectangle,  a  triangle,  and  a  trapezoid.  The  code  in 
subroutine  RESPNS  is  written  so  other  functions  may  be  easily  added,  and  the 
method  for  constructing  instrument  response  functions  is  described  in  the  source 
code.  When  a  waveband  calculation  is  performed  for  the  FGLOW  scenario,  the 
aerosol  optical  properties  used  are  for  the  center  wavelength  of  the  band.  If  the 
FGLOW  output  data  are  specified  as  apparent  temperature,  a  blackbody 
equivalent  temperature  is  calculated  based  on  a  conversion  of  the  band- 
integrated  radiation  to  a  weighted  average  value  at  the  center  wavelength  of  the 
band. 


55 


Specify  the  beginning  and  ending  wave  numbers  of  the  band  and  the  number 
of  calculation  points  in  the  band.  The  response  function  value  is  calculated  for 
each  of  the  points.  The  rectangle  has  values  of  1  in  the  band,  0  outside  the 
band,  and  0.5  at  the  boundaries.  The  triangle  has  a  peak  of  1  and  goes  to  0  at 
the  boundaries.  The  trapezoid  is  0  at  the  boundaries,  rises  linearly  from  0  to 
1  over  the  outer  20  percent  of  the  band  on  each  side,  and  has  a  value  of  1  for 
the  central  60  percent  of  the  band. 

The  number  of  calculation  points  in  the  band  must  be  odd  so  a  point  is  defined 
at  the  center  of  the  band.  The  code  checks  for  this  and  corrects  it  if  necessary. 
If  such  a  correction  is  made,  the  number  of  points  (and  therefore  the  wave 
number  increment)  differs  from  the  user  input. 


Table  11.  The  BAND  card 


Identifier  Variable 

BAND 

Default 

Description 

Wave  band  calculation  control 

IBAND 

0 

flag  wave  band  calculations: 

0  =»  single  wavelength  calculation 

1  =>  perform  wave  band  calculation 

WVNl 

beginning  wave  number  of  band  (cm  ') 

WVN2 

ending  wave  number  of  band  (cm  ') 

NWVN 

1 

number  of  wave  number  intervals  in  band 
(must  be  odd) 

IRESP 

0 

instrument  response  function  over  wave 
band: 

0  -  rectangular 

1  -  triangular 

2  -  trapezoidal 

IPSPEC 

0 

switch  printing  of  spectrally  resolved  values: 

0  =>  do  not  print  resolved  values 

1  ^  print  resolved  transmittance  and  radiance 
values 

56 


The  OPTl  card  permits  the  user  to  override  the  default  values  of  TCRITA,  the 
fractional  change  in  temperature  that  triggers  calls  to  ATLESF ;  NUMINT,  the 
number  of  calculation  intervals  through  the  plume  (should  be  an  even  number); 


and  NUMSEG,  the  number  of  calculation  intervals  between  the  plume  and 
observer.  The  default  values  speed  calculations  and  give  good  results. 

Table  12.  The  OPTl  card 

Identifier  Variable  Default  Description 

OPTl  Optional  parameters  for  calculation  control 


TCRITA 

NUMINT 

NUMSEG 


0.10  fractional  change  in  temperature  which 
triggers  a  call  to  the  band  model 
30  number  of  calculation  steps  through  the  plume 
1  number  of  path  segments  for  calculation  from 
plume  to  observer 


The  SCN3  card  supplies  imager  parameters  needed  for  scenario  3  calculations 
and  control.  Although  the  option  to  produce  apparent  temperature  output  image 
data  files  has  not  been  removed  from  the  code,  this  option  should  not  be  used. 
Accessory  programs,  described  in  appendix  B,  can  transform  the  radiance  image 
data  to  apparent  temperature  image  data. 

Table  13.  The  SCN3  card 

Identifier  Variable  Default  Description 

SCN3  Imager  parameters  for  scenario  3  (FGLOW) 


AINC 

resolution  per  pixel  (mrad) 

NGX 

number  of  columns  of  pixels 

NGZ 

number  of  rows  of  pixels 

lAPT 

0 

code  for  type  of  data  stored: 

0  =»  radiance  and  transmittance  data 

1  =►  apparent  temperature  data 

Use  the  TCAL  card  to  specify  an  internal  time-series  calculation.  If  FITTE  is 
called  by  other  modules  within  EOSAEL,  only  the  output  parameters  from  the 
final  time  are  passed  back  to  the  calling  routine.  Time-series  calculations 
suppress  printing  of  waveband  resolved  output.  All  times  are  converted  to 
integers  because  the  built-in  time  constant  for  calculations  is  1  s.  The  VELO 
parameter,  normally,  is  zero,  but  it  can  be  used  to  specify  movement  of  the 
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observer  toward  the  target  along  the  LOS  during  the  time  series  calculation. 
If  the  observer  crosses  the  model  system  wind-aligned  axis,  the  calculations  are 
aborted  and  a  message  is  written  to  the  output  file. 


Table  14.  The  TCAL  card 


Identifier  Variable 

TCAL 

Default 

Description 

Time-series  calculation  control 

STIME 

100 

start  time  for  calculations  (seconds  after 
start  of  fire) 

TINC 

1 

time  increment  between  calculations  (s) 

NCALC 

1 

number  of  calculations  in  time  series 

VELO 

0 

velocity  of  observer  toward  target  along 
LOS  (m/s) 

The  SVAR  card  lets  the  user  modify  certain  fire  parameters.  It  is,  primarily, 
intended  to  permit  matching  of  parameters  for  comparison  with  experimental 
data.  USE  THE  SVAR  CARD  WITH  CAUTION  —  IT  CAN  CAUSE 
UNREALISTIC  COMBINATIONS  OF  FIRE  PARAMETERS  AND 
ERRONEOUS  PREDICTIONS.  You  should  monitor  the  value  of  vertical  wind 
Wg  printed  in  the  output  file  when  the  fire  parameters  are  varied.  A  value  much 
outside  the  range  1.5  to  3  m/s  indicates  that  the  combined  fire  parameters  have 
changed  significantly  from  measured  values  and  the  model  predictions  may  not 
be  accurate.  See  appendix  C  for  a  discussion  of  varying  the  fire  parameters  and 
developing  model  fire  parameters  from  data. 


Table  15.  The  SVAR  card 


Identifier  Variable 

SVAR 

Default 

Description 

Variation  of  fire  parameters 

TEMPIN 

973.16 

fire  mean  temperature  (K) 

EFFAC 

1 

fractional  multiplier  of  aerosol  efficiency  factor 

RADIN 

radius  of  fire  (m) 

BTIME 

1500 

bum  time  (s) 
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The  mean  fire  temperature  TEMPIN  can  be  varied  fairly  fireely.  The  current 
value  seems  to  give  accurate  radiation  temperature  distributions  for  fires 
enclosed  within  metallic  hulls.  For  open  burning  of  fuel,  the  appropriate  mean 
temperature  may  be  up  to  several  hundred  degrees  higher  than  the  default  value 
used  in  the  model.  If  the  variation  of  the  fire  parameters  is  performed  to  match 
experimental  data,  the  peak  temperature  in  the  Gaussian  plume  at  the  fire  is 
roughly  2.9  times  the  mean  temperature. 

EFFAC  is  a  fi-actional  multiplier  that  allows  you  to  vary  the  creation  of  smoke 
aerosols.  Normally,  it  is  set  to  1.0,  but  can  be  any  number.  Values  between 
about  0.67  and  1.33  would  cover  the  general  range  of  smoke  variation  for 
hydrocarbon  fiiel  fires. 

RADIN  allows  you  to  vary  the  fire  radius.  Plume  scaling  in  the  atmosphere  is 
generally  considered  to  occur  for  fire  diameters  of  no  less  than  1  m,  so  the 
minimum  value  for  RADIN  should  be  0.5  m.  Note  that  the  actual  fire  image 
appears  to  have  approximately  two-thirds  the  specified  source  radius. 

BTIME  allows  you  to  vary  the  total  bum  time  for  the  fire.  The  current  value 
is  1500  s  (25  min).  Large  variations  fi'om  this  value  are  inappropriate  for  open 
burning  conditions,  but  substantial  increases  might  be  called  for  to  simulate  an 
enclosed  fire  in  which  the  bum  rate  is  controlled  by  an  inadequate  oxygen 
supply. 


Use  the  GO  card  to  indicate  that  repetitive  calculations  are  to  be  performed  for 
which  new  parameter  values  are  read  between  model  mns. 

Table  16.  The  GO  card 

Identifier  Variable  Default  Description 

GO  Execute  and  read  data  for  a  subsequent  mn 
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Use  the  DONE  card  to  execute  a  single  run  or  the  last  of  a  series  of  runs.  It 
returns  control  to  the  driver  program. 

Table  17.  The  DONE  card 

Identifier  Variable  Default  Description 

DONE  Execute  and  return  control  to  EOEXEC 


The  default  values  shown  on  the  cards  are  used  if  no  value  is  furnished.  As  can 
be  seen,  not  all  variables  have  default  values.  Other  (default  substitute)  values 
are  used  when  certain  input  parameters  are  outside  an  appropriate  range  of 
values.  A  notice  that  identifies  the  input  card  involved  is  printed  in  the  output 
file  if  any  substitution  of  values  is  made. 

Selecting  a  turbulence  effects  calculation,  a  time-series  calculation,  or  a 
scenario  3  (FGLOW)  calculation  (particularly  for  waveband  calculations)  may 
increase  run  time  significantly.  Typical  run  times  on  a  20-MHz  386  PC  with 
coprocessor  are  shown  in  table  18. 


Table  18.  Run  times  for  selected  run  conditions 


single  LOS 

single  wavelength 

1  min 

12  s 

single  LOS 

25  wave  number  band 

0  min 

13  s 

single  LOS 

single  wavelength 

25  point  time  series 

with  turbulence  effects 

14  min 

23  s 

without  turbulence  effects 

0  min 

34  s 

single  LOS 

25  wave  number  band 

1  min 

2  s 

25  point  time  series 

25  by  25  array 

single  wavelength 

3  min 

45  s 

(625  total  LOS) 

25  by  25  array 

25  wave  number  band 

39  min 

28  s 

(625  total  LOS) 
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4.2  Output 


4.2.1  Parameters  Returned  to  the  Calling  Routine 

The  parameters  available  for  return  to  the  calling  routine  were  increased.  They 
now  include  the  following: 

FITRAN  -  the  total  transmittance  for  the  LOS 
lERR  -  a  code  to  indicate  an  invalid  input  card  identifier 
ATRAD  -  the  attenuated  target  radiance  for  the  LOS 
PATRAD  -  the  emitted  radiance  firom  the  plume  for  the  LOS 
CONLEN  -  the  aerosol  concentration-length  product  for  the  LOS 
SPTMAG  -  the  turbulence-induced  spot  (area)  magnification 
JITRAD  -  the  turbulence-induced  beam  wander  radius. 

All  the  parameters  (except  the  error  code)  are  output  for  scenario  1  or  2 
calculations  that  include  turbulence  effects.  If  turbulence  effects  calculations 
are  omitted,  the  last  two  parameter  values  are  meaningless. 

For  FGLOW  calculations,  the  output  parameters  are  those  for  the  LOS  with  the 
highest  emitted  radiance. 

4.2.2  Output  Printed  or  Stored  in  Files 

FITTE  scenario  1  or  2  output  is  contained  in  a  single  file.  Samples  are  shown 
in  the  example  section.  Most  of  the  input  data  are  listed,  plume  parameters  are 
given  for  the  calculation  point  on  the  LOS  closest  to  the  plume  centerline,  and 
the  model  predictions  are  printed.  The  output  parameters  provided  are  modeled 
calculation  time,  total  transmittance,  aerosol  transmittance,  attenuated  target 
radiance,  emitted  radiance,  apparent  temperature,  concentration-length,  and,  for 
single  wavelength  calculations  that  include  turbulence  effects,  laser  spot 
magnification,  and  jitter  radius. 

If  transmittance  and  radiance  are  calculated  over  a  spectral  band,  you  have  the 
option  (set  IPSPEC  =  1  on  BAND)  of  printing  the  spectrally  resolved  values 
of  the  radiative  transfer  predictions.  Quantities  printed  are  the  wavelength, 
wave  number,  transmittance,  target  radiance,  emitted  radiance,  total  radiance. 
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and  instrument  spectral  response.  You  cannot  exercise  this  option  if  a  time- 
series  or  an  FGLOW  calculation  is  being  performed. 

For  the  FGLOW  scenario,  two  output  files  are  created.  The  standard  output  file 
contains  a  summary  of  input  data  and  information  about  the  matrix  data  file 
created.  The  second  file  created  contains  a  header  with  information  on 
calculation  parameters  and  ngx  times  ngz  rows  that  contain  emitted  radiance  and 
total  transmittance  data  or  apparent  temperature  data  for  the  image  matrix 
specified  on  the  SCN3  card.  The  examples  show  samples  for  running  FITTE 
for  various  conditions.  The  output  file  is  named  FGLOW l.FDT  for  a  single 
image  calculation.  If  a  time-series  of  n  images  is  calculated,  the  files  are 
named  FGLOWl.FDT  to  FGLOWn.FDT.  The  FGLOW  data  files  should  be 
renamed  if  you  wish  to  keep  them  because  the  model  overwrites  one  (or  more) 
of  the  files  the  next  time  it  is  run. 


5.  Sample  Runs 

This  section  shows  examples  for  running  FITTE  for  various  conditions.  The 
examples  show  the  input  card  and  output  files  and  may  be  accompanied  by 
comments,  an  FGLOW  data  file,  and  illustrations  of  FGLOW  output. 

The  examples  and  output  were  current  at  the  time  the  manual  was  prepared  for 
publication.  Some  changes  may  occur  because  of  code  modifications.  Check 
the  sample  input/ouQiut  files  from  the  EOSAEL  tape  to  identify  any  such 
changes. 

The  examples  were  chosen  to  illustrate  most  of  the  model  options  and  were  run 
under  the  following  conditions: 

(1)  The  synthetic  wind  fluctuation  file  FITEX.DAT  was  copied  to 
FIT4D.DAT.  To  obtain  more  realistic  output,  one  of  the  wind  fluctuation  files 
based  on  field  data  (see  information  in  the  file  README.FIT)  should  be  copied 
to  FIT4D.DAT  after  the  examples  are  run. 

(2)  The  input  files  listed  are  for  use  with  the  PC  version  of  the  code  using 
FITDRV  as  the  main  program.  Several  additional  cards  are  needed  to  use  the 
EOSAEL  main  driver.  The  additional  cards  include  a  WAVL  and  a  FITTE  card 
as  a  minimum,  and  may  include  cards  with  climate  and  geometry  data. 

(3)  The  EOSAEL  header  page  was  omitted  from  the  output  files. 

(4)  The  DONE,  END,  and  STOP  cards  were  omitted  from  most  of  the 
examples. 
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5.1  Example  1:  A  Scenario  1  Calculation  for  a  Single 
Wavelength 

The  user  coordinate  system  positive  x-axis  points  east,  and  a  type  2  fire  is  used 
for  a  scenario  1  calculation  through  a  4D  plume.  The  target  temperature  is 
10  °C  and  has  an  emissivity  of  0.8.  The  initial  laser  beam  diameter  is  0.3  cm. 


INPUT: 

REFD 

90.0 

2. 

1. 

1. 

SCEN 

DETD 

8.0 

0.30 

-100.0 

6.0 

8.0 

900.0 

6.0 

METD 

SRCD 

3.0 

0.0 

270. 

0.0 

6.5 

0.0 

83.0 

1200.0 

3.0 

TARG  10.  0.8 

DONE 

END 

STOP 

OUTPUT: 

★4r**4r**4r***4r****** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 

★***4r4r**4r4r**iilr4r**** 


FIRE  TYPE:  TRUCK 


X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  M 

( 

X  , 

Y 

z 

) 

OBSERVER 

( 

8.0, 

-100.0, 

6.0 

) 

TARGET 

( 

8.0, 

900.0, 

6.0 

) 

FIRE 

< 

0.0, 

0.0, 

0.0 

) 

METEOROLOGICAL  DATA  : 


WINDSPEED 
TEMPERATURE 
AIR  DENSITY 

3.0  M/S 

6.5  C 

1200.0  G/M**3 

WIND  DIRECTION 
RELATIVE  HUMIDITY 
PASQUILL  CATEGORY 

270.0  DEG 
83.0  % 

3 

TARGET  TEMPERATURE  10.0  C 

EMISSIVITY 

0.80 

LASER  WAVELENGTH  0.63  UM 

BEAM  DIAMETER 

0.30  CM 

PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.39  M/S  COEFFICIENT  (C)  6.906 
BUOYANCY  FLUX  (F)  37.5  VIRTUAL  SOURCE  POSITION  (M) 
ENTRAINMENT  CONSTANT  0.4134  XV:  -7.12,  ZV:  -8.52 


PLUHE-LOS  GEOMETRIC  AND  DISTRIBUTION  FUNCTION  DATA  AT  CLOSEST  APPROACH: 

CENTERLINE  COORDINATES  (  7.26,  -0.12,  5.68  )  M 

LOS  COORDINATES  (  8.00,  -0.24,  6.00  )  M 

DISTANCE  0.82  M  EXCESS  TEMPERATURE  153.  C 

PLUME  RADIUS  2.32  M  AEROSOL  DENSITY  0.25  G/H**3 

****************************** 

*  4D  FITTE  PREDICTIONS  * 
****************************** 

TRANSMITTANCE  ATTENUATED  SPOT  JITTER 

TIME  TOTAL  AEROSOL  TARG  RAD  PATH  RAD  APP  TEMP  CL  HAGN  RADIUS 
S6C  watts/ni**2/sr/uni  kelvin  g/m**2  m 

100.  0.00  0.00  0.129E-28  0.275E-15  403.0  1.25  11.6  0.00 
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5.2 


Example  2:  A  Scenario  2  Calculation  for  a  Single 
Wavelength 


The  only  change  from  example  1  is  that  the  calculation  is  done  for  scenario  2. 
Note  the  effect  of  reversing  the  direction  of  the  calculation  on  the  turbulence 
effects  output. 


REFD 

90,0 

2. 

2. 

1. 

SCEN 

8.0 

-100.0 

6.0 

8.0 

900.0 

6.0 

DETD 

0.30 

METD 

3.0 

270. 

6.5 

83.0 

1200.0 

3.0 

SRCD 

0.0 

0.0 

0.0 

TARG 

10. 

0.8 

OUTPUT: 

****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X-AXIS  HEADING:  90. 

0  DEG  CW  FROM  NORTH 

POSITION  M  (  X  , 

Y  ,  2  ) 

OBSERVER  (  8.0, 

-100.0,  6.0  ) 

TARGET  (  8.0, 

900.0,  6.0  ) 

FIRE  (  0.0, 

0.0,  0.0  ) 

METEOROLOGICAL  DATA  : 

WINDSPEED  3.0  M/S 

WIND  DIRECTION 

270.0  DEG 

TEMPERATURE  6.5  C 

RELATIVE  HUMIDITY 

83.0  % 

AIR  DENSITY  1200.0  G/M**3 

PASQUILL  CATEGORY 

3 

TARGET  TEMPERATURE  10.0  C 

EMISSIVITY  0.80 

LASER  WAVELENGTH  0.63  UM 

BEAM  DIAMETER 

0.30  CM 

PLUME  CONSTANTS  : 

VERTICAL 

WIND 

(WO) 

2.39  M/S 

COEFFICIENT  (C)  6.906 

BUOYANCY 

FLUX 

(F) 

37.5 

VIRTUAL  SOURCE  POSITION  (M) 

ENTRAINMENT  CONSTANT 

0.4134 

XV:  -7.12,  ZV:  -8.52 

PLUME-LOS  GECWETRIC  AND  DISTRIBUTION  FUNCTION  DATA  AT  CLOSEST  APPROACH: 


CENTERLINE  COORDINATES  <  7.26,  -0.12, 
LOS  COORDINATES  (  8.00,  -0.24, 
DISTANCE  0.82  M  EXCESS  TEMPERATURE 
PLUME  RADIUS  2.32  M  AEROSOL  DENSITY 

5.68  )  M 
6.00  )  M 
153.  C 
0.25  G/M**3 

****************************** 

*  4D  FITTE  PREDICTIONS  * 

****************************** 

TIME 

sec 

TRANSMITTANCE 
TOTAL  AEROSOL 

ATTENUATED 

TARG  RAD  PATH  RAD  APP  TEMP 

watts/m**2/sr/um  kelvin 

CL 

g/m**2 

SPOT 

MAGN 

JITTER 

RADIUS 

m 

100. 

0.00  0.00 

0.129E-28  0.275E-15  403.0 

1.25 

1.0 

0.00 
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5.3  Example  3:  A  Scenario  1  Calculation  for  a  (Visible) 
Waveband 

FITTE  does  not  include  molecular  effects  for  this  waveband. 


INPUT:  REFD 

90.0 

2. 

1. 

1. 

SCEN 

8.0 

-100.0 

6.0 

8.0 

900.0 

6.0 

DETD 

0.30 

METD 

3.0 

270. 

6.5 

83.0 

1200.0 

3.0 

SRCD 

0.0 

0.0 

0.0 

TARG 

10. 

0.8 

BAND 

1. 

13000. 

19231. 

25. 

2. 

0. 

OUTPUT : 

****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X“AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  M 

( 

X  , 

Y 

2 

)  : 

OBSERVER 

( 

8.0, 

-100.0, 

6.0 

) 

TARGET 

( 

8.0, 

900.0, 

6.0 

) 

FIRE 

( 

0.0, 

0.0, 

0.0 

) 

WINDSPEED 
TEMPERATURE 
AIR  DENSITY 


METEOROLOGICAL  DATA  : 

3.0  M/S  WIND  DIRECTION 

6.5  C  RELATIVE  HUMIDITY 

1200.0  G/M**3  PASQUILL  CATEGORY 


270.0  DEG 
83.0  % 

3 


TARGET  TEMPERATURE  10.0 


EMISSIVITY  0.80 


WAVEBAND  DATA  : 

WAVEBAND  FROM  13000.00  CM-1  TO  19231.00  CM-1  CENTER  WAVELENGTH  0.621  UM 
25  POINTS  RESOLUTION  INSTRUMENT  RESPONSE:  TRAP  MOLECULES  INCLUDED:  YES 
WAVEBAND  DOES  NOT  OVERLAP  FITTE  MOLECULAR  DATA 


VERTICAL  WIND  (WO) 
BUOYANCY  FLUX  (F) 
ENTRAINMENT  CONSTANT 


PLUME  CONSTANTS  : 

2.39  M/S  COEFFICIENT  (C)  6.906 
37.5  VIRTUAL  SOURCE  POSITION  (M) 
0.4134  XV:  -7.12,  ZV:  -8.52 


PLUME-LOS  GEOMETRIC  AND  DISTRIBUTION  FUNCTION  DATA  AT  CLOSEST  APPROACH: 


CENTERLINE  COORDINATES 
LOS  COORDINATES 
DISTANCE  0.82  M 

PLUME  RADIUS  2.32  M 


(  7.26,  -0.12,  5.68  )  M 

(  8.00,  -0.24,  6.00  )  M 

EXCESS  TEMPERATURE  153.  C 

AEROSOL  DENSITY  0.25  G/M**3 


****************************** 

*  4D  FITTE  PREDICTIONS  * 
****************************** 


TRANSMITTANCE  ATTENUATED 


TIME 

TOTAL 

AEROSOL 

TARG  RAD 

PATH  RAD 

APP  TEMP 

CL 

sec 

watts/m**2/sr 

kelvin 

g/m**2 

100. 

0.00 

0.00 

0.283E-36 

0.860E-14 

448.5 

1.25 
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5.4 


Example  4:  A  Scenario  1  Time-Series  Calculation  for  a 
Single  Wavelength 


The  TCAL  card  specifies  a  time-series  calculation.  The  calculation  starts  100  s 
after  the  start  of  the  fire.  The  interval  between  calculations  is  1  s,  and  1 1 
calculations  are  performed  to  span  a  10-s  period.  At  this  wavelength 
(1.06  pm),  the  inclusion  of  molecular  absorption  has  no  effect  on  the 
transmittance. 


INPUT: 


REFD 

90.000 

3.000 

1. 

1. 

SCEN 

5.000 

400.000 

3.00 

5.000 

-800.000 

3.0 

SRCD 

0.000 

0.000 

0.000 

DETD 

5.000 

TARG 

37.0 

1.0 

METD 

4.000 

270.000 

27.000 

60.000 

1009.728 

4.000 

TCAL 

100. 

1. 

11. 

0. 

MOLS 

1. 

1. 

CXJTPUT: 

****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TANK 


X-AXIS  HEADING: 

90.0  DEG  CU  FROM 

NORTH 

POSITION  M 

( 

X  , 

Y 

Z 

) 

OBSERVER 

( 

5.0, 

400.0, 

3.0 

) 

TARGET 

( 

5.0, 

-800.0, 

3.0 

) 

FIRE 

( 

0.0, 

0.0, 

0.0 

) 

UINDSPEED 
TEMPERATURE 
AIR  DENSITY 


METEOROLOGICAL  DATA  : 

4.0  M/S  WIND  DIRECTION 


27-0  C 

1009.7  G/M**3 


RELATIVE  HUMIDITY 
PASQUILL  CATEGORY 


270.0  DEG 
60.0  % 

4 


TARGET  TEMPERATURE  37.0 


EMISSIVITY  1.00 


LASER  WAVELENGTH  1.06  UM  BEAM  DIAMETER  5.00  CM 


VERTICAL  WIND  (WO) 
BUOYANCY  FLUX  (F> 
ENTRAINMENT  CONSTANT 


PLUME  CONSTANTS  : 

2.78  M/S  COEFFICIENT  (C)  10.887 

169.1  VIRTUAL  SOURCE  POSITION  (M) 

0.4433  XV:  -17.71,  ZV:  -18.49 


PLUME-LOS  GEOMETRIC  AND  DISTRIBUTION  FUNCTION  DATA  AT  CLOSEST  APPROACH: 


CENTERLINE  COORDINATES 
LOS  COORDINATES 
DISTANCE  1.15  M 

PLUME  RADIUS  2.87  M 


(  3.87,  -0.52, 

(  5-00,  -0.50, 

EXCESS  TEMPERATURE 
AEROSOL  DENSITY 


3.24  )  M 
3.00  )  M 
0.  C 
1.28  G/M**3 
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itifkititifk*1t******iticit*it********** 

*  4D  FITTE  PREDICTIONS  * 


FITTE  OUTPUT  FOR  TIME-SERIES  CALCULATION 


TRANSMITTANCE 

ATTENUATED 

SPOT 

JITTER 

TIME 

TOTAL 

AEROSOL 

TARG  RAD 

PATH  RAD 

APP  TEMP 

CL 

MAGN 

RADIUS 

sec 

watts/m**2/sr/uni 

keLvin 

g/m**2 

m 

100. 

0.00 

0.00 

0.107E-14 

0.296E-01 

622.6 

3.05 

172.4 

0.44 

101. 

0.00 

0.00 

0.446E-14 

0.321E+01 

792.9 

2.56 

224.9 

0.49 

102. 

0.00 

0.00 

0.317E‘17 

0.157E+00 

674.2 

5.03 

477.0 

0.66 

103. 

0.00 

0.00 

0.278E-13 

0.134E-03 

499.1 

1.94 

80.6 

0.32 

104. 

0.00 

0.00 

0.250E-15 

0.304E+00 

697.0 

3.54 

539.4 

0.69 

105. 

0.00 

0.00 

0.832E-15 

0.596E+00 

721.9 

3.14 

327.5 

0.56 

106. 

0.00 

0.00 

0.171E-15 

0.707E-01 

648.5 

3.67 

291.5 

0.54 

107. 

0.05 

0.05 

0.409E-12 

0.148E’04 

461.8 

1.03 

49.6 

0.26 

108. 

0.06 

0.06 

0.530E-12 

0.159E-04 

462.9 

0.94 

61.6 

0.29 

109. 

0.03 

0.03 

0.240E-12 

0.688E-03 

531.1 

1.21 

70.9 

0.30 

110. 

0.00 

0.00 

0.727E-14 

0.865E+00 

736.5 

2.40 

340.0 

0.57 
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5.5  Example  5:  A  Scenario  1  Calculation  With  Variation  of  Fire 
Parameters  and  Printout  of  Spectrally  Resolved  Results  for 
a  Selected  Instrument  Response  Function 

This  calculation  is  performed  over  an  IR  waveband.  The  last  two  parameters 
on  the  BAND  card  specify  use  of  a  trapezoidal  response  function  and  printout 
of  spectrally  resolved  output.  Note  the  use  of  the  SVAR  card  to  specify  a 
higher  fire  temperature. 


90.0  2.  1.  1. 

0.0  0.0  80.0  30.0  1800.0  20.0 

0.30 

3.0  255.  6.5  83.0  1200.0  3.0 

-6.0  1650.0  0.0 

10.  0.8 

1223.  1. 

1.  2600.  3300.  25.  2.  1. 


****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 
USER  ENTERED  FIRE  VARIATION: 

AVERAGE  TEMPERATURE:  1223.  K 

X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 

POSITION  M  (  X  ,  Y  ,  Z  )  : 

OBSERVER  (  0.0,  0.0,  80.0  ) 

TARGET  (  30.0,  1800.0,  20.0  ) 

FIRE  (  -6.0,  1650.0,  0.0  ) 

METEOROLOGICAL  DATA  : 

WINDSPEED  3.0  M/S  WIND  DIRECTION  255.0  DEG 

TEMPERATURE  6.5  C  RELATIVE  HUMIDITY  83.0  % 

AIR  DENSITY  1200.0  G/M**3  PASQUILL  CATEGORY  3 

TARGET  TEMPERATURE  10.0  C  EMISSIVITY  0.80 

WAVEBAND  DATA  : 

WAVEBAND  FROM  2600.00  CM-1  TO  3300.00  CM-1  CENTER  WAVELENGTH  3.390  UM 
25  POINTS  RESOLUTION  INSTRUMENT  RESPONSE:  TRAP  MOLECULES  INCLUDED:  YES 

PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.16  M/S  COEFFICIENT  (C)  6.666 

BUOYANCY  FLUX  (F)  37.5  VIRTUAL  SOURCE  POSITION  (M) 

ENTRAINMENT  CONSTANT  0.4359  XV:  -8.68,  ZV:  -9.38 


PLUME-LOS  GEOMETRIC  AND  DISTRIBUTION  FUNCTION  DATA  AT  CLOSEST  APPROACH: 

CENTERLINE  COORDINATES  <  25.81,  1657.55,  19.09  )  H 
LOS  COORDINATES  (  27.65,  1659.16,  24.69  )  M 
DISTANCE  6.12  M  EXCESS  TEMPERATURE  0.  C 
PLUME  RADIUS  7.24  M  AEROSOL  DENSITY  0.00  G/M**3 


INPUT:  REFD 

SCEN 
DETD 
METD 
SRCD 
TARG 
SVAR 
BAND 

OUTPUT: 
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****************************** 

*  4D  FITTE  PREDICTIONS  * 
****************************** 


TRANSMITTANCE 
TIME  TOTAL  AEROSOL 
sec 

100.  0.85  0.85 


ATTENUATED 
TARG  RAD  PATH  RAD 
watts/m**2/sr 

0.105E-01  0.714E-02 


APP  TEMP  CL 
kelvin  g/m**2 

263.9  0.12 


SPECTRAL  RADIANCE  AND  TRANSMITTANCE 


WAVELTH 

FREQ 

(UM) 

(CM**-1) 

3.8462 

2600.000 

3.8035 

2629.167 

3.7618 

2658.333 

3.7209 

2687.500 

3.6810 

2716.667 

3.6419 

2745.833 

3.6036 

2775.000 

3.5661 

2804.167 

3.5294 

2833.333 

3.4934 

2862.500 

3.4582 

2891.667 

3.4237 

2920.833 

3.3898 

2950.000 

3.3566 

2979.167 

3.3241 

3008.333 

3.2922 

3037.500 

3.2609 

3066.667 

3.2301 

3095.833 

3.2000 

3125.000 

3.1704 

3154.167 

3.1414 

3183.333 

3.1128 

3212.500 

3.0848 

3241.667 

3.0573 

3270.833 

3.0303 

3300.000 

TRANS 

0.8556 

0.8552 

0.8544 

0.8535 

0.8527 

0.8519 

0.8511 

0.8503 

0.8496 

0.8488 

0.8481 

0.8474 

0.8467 

0.8460 

0.8453 

0.8447 

0.8440 

0.8434 

0.8428 

0.8422 

0.8416 

0.8410 

0.8404 

0.8399 

0.8393 


TARG  RAD  PATH  RAD  TOTAL  RAD  INSTR 
(WATTS/M**2/SR/CM)  RESPONSE 


2.423E-05  4.050E-05 
2.369E-05  3.616E-05 
2.315E-05  3.237E-05 
2.262E-05  2.897E-05 
2.212E-05  2.592E-05 
2.163E-05  2.318E-05 
2.116E-05  2.071E-05 
2.070E-05  1.850E-05 
2.026E-05  1.652E‘05 
1.983E-05  1.475E-05 
1.942E-05  1.316E-05 
1.902E-05  1.173E-05 
1.863E-05  1.046E-05 
1.825E-05  9.319E-06 
1.788E-05  8.300E’06 
1.753E-05  7.390E“06 
1.718E-05  6.578E-06 
1.685E-05  5.852E-06 
1.652E-05  5.205E-06 
1.621E-05  4.628E-06 
1.590E-05  4.113E-06 
1.560E-05  3.655E-06 
1.531E-05  3.246E-06 
1.503E-05  2.883E-06 
1.476E-05  2.559E’06 


6.473E-05 

0.0000 

5.984E-05 

0.2083 

5.552E-05 

0.4167 

5.160E-05 

0.6250 

4.804E-05 

0.8333 

4.481E-05 

1.0000 

4.187E-05 

1.0000 

3.921E-05 

1.0000 

3.678E-05 

1.0000 

3.458E-05 

1.0000 

3.257E‘05 

1.0000 

3.075E-05 

1.0000 

2.908E-05 

1.0000 

2.757E-05 

1.0000 

2.618E-05 

1.0000 

2.492E-05 

1.0000 

2.376E-05 

1.0000 

2.270E-05 

1.0000 

2.173E-05 

1.0000 

2.083E-05 

1.0000 

2.001E-05 

0.8333 

1.926E-05 

0.6250 

1.856E-05 

0.4167 

1.791E-05 

0.2083 

1.731E’05 

0,0000 
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5.6  Example  6:  A  Single  Wavelength  Scenario  3  Calculation 


The  SCN3  card  provides  all  parameters  except  wavelength  and  waveband  for 
the  scenario  3  imager.  The  angular  resolution  is  0.5  mrad  per  pixel.  A 
5-column,  3 -row  array  of  radiance  and  temperature  is  stored  after  the  parameter 
summary  line  in  the  data  file  (shown  on  the  next  page).  The  pixel  size  is  large 
enough,  compared  with  the  plume  width,  that  the  values  in  the  data  array  are 
area  averaged.  The  -1  on  the  REFD  card  specifies  a  mean-value  plume.  The 
FGLOWl.FDT  file  is  shown  on  the  next  page. 


INPUT: 

REFD 

SCEN 

90.000 

25.000 

.  2.000 
-1000.000 

3. 

20.00 

-1. 

25.000 

100.000 

20.0 

SRCD 

METD 

0.000 

3.000 

0.000 

270.000 

0.000 

27.000 

60.000 

1009.728 

4.000 

SCN3 

0.50 

5. 

3. 

0. 

OUTPUT: 

****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X-AXIS  HEADING;  90.0  DEG  CU  FROM  NORTH 


POSITION  M  ( 
OBSERVER  ( 
TARGET  ( 
FIRE  ( 


X  ,  Y  ,  Z  ) 
25.0,  -1000.0,  20.0  ) 

25.0,  100.0,  20.0  ) 

0.0,  0.0,  0.0  ) 


WINDSPEED 
TEMPERATURE 
AIR  DENSITY 


METEOROLOGICAL  DATA  : 

3.0  M/S  WIND  DIRECTION  270.0  DEG 

27.0  C  RELATIVE  HUMIDITY  60.0  % 

1009.7  G/H**3  PASQUILL  CATEGORY  4 


IMAGER  DATA: 

TOTAL  FOV:  2.50  BY  1.50  (MR)  RESOLUTION  PER  PIXEL  0.50  MR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  5,  3  ) 

WAVELENGTH  10.60  UM 


VERTICAL  WIND  (WO) 
BUOYANCY  FLUX  (F) 
ENTRAINMENT  CONSTANT 


PLUME  CONSTANTS  : 

2.74  M/S  COEFFICIENT  (C)  7.530 

41.5  VIRTUAL  SOURCE  POSITION  (M) 

0.3820  XV:  -6.17,  ZV:  -8.45 


**************** 

*  FGLOW  OUTPUT  * 
**************** 


AT  100.  SECONDS  AFTER  START  OF  FIRE 


RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR 
THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  24.00, 

AND  THE  ARRAY  CENTER  IS  AT  (  25.00, 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT 


5  COLUMN,  3  ROW  ARRAY. 

0.00,  20.50  )  M 

0.00,  20.00  )  M 

THE  SOURCE  DISTANCE 


FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.500  M  UNITS:  Watts/in**2/sr/um 


VALUES  ARE  AREA  AVERAGED 
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The  data  file  is  listed  below.  The  first  row  contains  the  following  information: 

•  the  (single  or  center-of-waveband)  wavelength  for  the  calculation 

•  the  number  of  columns  and  rows  of  imager  data 

•  the  coordinates  of  the  center  of  the  upper  left  pixel  in  the  plane  of  the  fire 

•  the  angular  resolution  for  one  pixel 

•  the  coordinates  of  the  center  pixel  in  the  array 

•  the  ambient  temperature 

•  a  waveband  function  (used  to  convert  waveband  integrated  radiance  to  a 
single- wavelength  blackbody  equivalent  temperature). 


The  succeeding  rows  each  contain  (radiance  and  transmittance)  data  for  one 
pixel. 


10.60000  5 

2.071363E+00 
2.438133E+00 
2.823721E+00 
3.216068E+00 
3.610541E+00 
2.705834E+00 
3.116120E+00 
3.534536E+00 
3.946305E+00 
4.343624E+00 
3.436608E+00 
3.874342E+00 
4.302427E+00 
4.709477E+00 
5.090297E+00 


3  24.000 

0.774 
0.738 
0.701 
0.666 
0.631 
0.713 
0.675 
0.638 
0.603 
0.571 
0.648 
0.611 
0.576 
0.544 
0.515 


0.000  20.500 


0.500  25.000 


0.000  20.000  300.16  1.00 
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5.7  Example  7:  A  Scenario  3  Waveband  Calculation 

For  this  example,  the  calculation  is  done  for  a  mean-value  plume.  The  fourth 
parameter  on  the  SCN3  card  is  set  to  1  to  produce  apparent  temperature  rather 
than  radiance  and  transmittance  ou^ut.  The  FGLOWl.FDT  file  is  shown  on 
the  next  page.  See  example  6  for  a  description  of  the  information  included  in 
the  file. 

This  is  no  longer  a  recommended  mode  of  operation.  Unless  apparent 
temperature  information  is  the  only  information  desired,  it  is  more  efficient  to 
create  a  radiance  and  transmittance  file  and  process  it  with  the  accessory 
programs  described  in  appendix  B  to  obtain  apparent  temperature  output. 


REFD 

SCEN 

90.000 

25.000 

2.000 

-1000.000 

3. 

20.00 

-1. 

25.000 

100.000 

20.0 

SRCD 

METD 

0.000 

3.000 

0.000 

270.000 

0.000 

27.000 

60.000 

1009.728 

4.000 

SCN3 

BAND 

0.50 

1. 

5. 

793.65 

3. 

1162.79 

1. 

15. 

2. 

0. 

OUTPUT: 

****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 

X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  H  ( 
OBSERVER  ( 
TARGET  ( 
FIRE  ( 


X  ,  Y  ,  Z  ) 
25.0,  -1000.0,  20.0  ) 

25.0,  100.0,  20.0  ) 

0.0,  0.0,  0.0  ) 


WINDSPEED 
TEMPERATURE 
AIR  DENSITY 


METEOROLOGICAL  DATA  : 

3.0  M/S  WIND  DIRECTION  270.0  DEG 

27.0  C  RELATIVE  HUMIDITY  60.0  % 

1009.7  G/M**3  PASQUILL  CATEGORY  4 


IMAGER  DATA: 

TOTAL  FOV:  2.50  BY  1.50  (MR)  RESOLUTION  PER  PIXEL  0.50  MR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  5,  3  ) 

WAVEBAND  FROM  793.65  CM-1  TO  1162.79  CM-1  CENTER  WAVELENGTH  10.223  UM 
15  POINTS  RESOLUTION  INSTRUMENT  RESPONSE:  TRAP  MOLECULES  INCLUDED:  YES 

APPARENT  TEMPERATURE  CALCULATED  FOR  CENTER  WAVELENGTH 


VERTICAL  WIND  (WO) 
BUOYANCY  FLUX  (F) 
ENTRAINMENT  CONSTANT 


PLUME  CONSTANTS  : 

2.74  M/S  COEFFICIENT  (C)  7.530 

41.5  VIRTUAL  SOURCE  POSITION  (M) 

0.3820  XV:  -6.17,  ZV:  -8.45 
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**************** 

*  FGLOW  OUTPUT  * 

*4r4r4r4r*********** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

APPARENT  TEMPERATURE  DATA  STORED  FOR  5  COLUMN,  3  ROW  ARRAY. 
THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  24.00,  0.00,  20.50  )  M 

AND  THE  ARRAY  CENTER  IS  AT  (  25.00,  0.00,  20.00  )  M 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SOURCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.500  M  UNITS:  degrees  K 

VALUES  ARE  AREA  AVERAGED 


The  FGLOW1.FDT  file  is  shown  be low - 

10.22265  5  3  24.000  0.000  20.500 

299. 

299. 

299. 

299. 

299. 

299. 

299. 

299. 

300. 

300. 

299. 

300. 

300. 

301. 

301. 


0.500  25.000  0.000  20.000  300.16  0.31 
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5.8  Example  8:  A  Time-Series  of  Single  Wavelength  Scenario  3 
Calculations 


INPUT: 


OUTPUT: 


The  SCN3  and  TCAL  cards  are  used  to  create  a  sequence  of  images.  The 
output  data  files,  FGLOWl  .FDT  through  FGLOW6.FDT,  are  not  shown  for  this 
example.  The  output  creates  a  series  of  files  separated  by  1  min  that  span  a 
5-min  period.  Figures  10a  and  b  show  contour  plots  of  the  data.  The  accessory 
program  XS2GRD  was  used  to  transform  the  data  files  to  input  files  for  the 
plotting  program. 


REFD 

90.0 

2. 

3. 

1. 

SCEN 

5.0 

-1000.0 

6.0 

5.0 

100.0 

6.0 

SRCD 

0.0 

0.0 

0.0 

METD 

3.0 

270.0 

27.0 

60.0 

1009.728 

4.0 

DETD 

0.125 

TARG 

10.0 

0.8 

SCN3 

0.125 

121. 

81. 

0. 

TCAL 

100. 

60. 

6. 

0. 

****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X-AXIS 

HEADING: 

90.0  DEG  CW 

FROM  NORTH 

POSITION  M  ( 

X  . 

Y 

,  z  ) 

OBSERVER  ( 

5.0, 

-1000.0, 

6.0  ) 

TARGET 

( 

5.0, 

100.0, 

6.0  ) 

FIRE 

( 

0.0, 

0.0, 

0.0  ) 

METEOROLOGICAL  DATA 

: 

WINDSPEED 

3.0  M/S 

WIND  DIRECTION 

TEMPERATURE 

27.0  C 

RELATIVE 

HUMIDITY 

AIR  DENSITY 

1009.7  G/M**3 

PASQUILL 

CATEGORY 

TOTAL  FOV: 


270.0  DEG 
60.0  % 

4 


IMAGER  DATA: 

15.13  BY  10.13  (MR)  RESOLUTION  PER  PIXEL  0.13  MR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  121,  81  ) 

WAVELENGTH  10.60  UM 


VERTICAL  WIND  (WO) 
BUOYANCY  FLUX  (F) 
ENTRAINMENT  CONSTANT 


PLUME  CONSTANTS  : 

2.74  M/S  COEFFICIENT  (C)  7.530 
41.5  VIRTUAL  SOURCE  POSITION  (M) 
0.3820  XV:  -6.17,  ZV:  -8.45 


**************** 

*  FGLOW  OUTPUT  * 

**************** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR  121  COLUMN,  81  ROW  ARRAY. 

THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  -2.50,  0.00,  11.00  )  M 

AND  THE  ARRAY  CENTER  IS  AT  (  5.00,  0.00,  6.00  )  M 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SOURCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.125  M  UNITS:  Watts/in**2/sr/utn 

THIS  WAS  REPEATED  6  TIMES  AT  60  SECOND  INTERVALS 
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Figure  10a.  Transmittance  plots  created  with  data  from  the  example 
8  simulation.  Wavelength  was  10.6  /tm.  Tic  spacing  is  5  m. 


Figure  10b.  Radiance  plots  created  with  data  from  the  example  8 
simulation.  Wavelength  was  10.6  /tm.  Tic  spacing  is  5  m. 
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5.9  Example  9:  Internal  Cycling  With  Variation  of  Wind 
Direction  to  Produce  Five  Output  Files 

Example  9  illustrates  the  use  of  the  GO  and  DONE  cards  to  perform  a  series 
of  calculations  with  variation  of  parameters  other  than  time.  The  five  pages 
of  the  EX9.0UT  file  are  shown.  Figure  11  was  prepared  fi-om  the  five 
(FGLOWn.FDT)  data  files.  The  accessory  program  XS2GRD  transformed  the 
data  files  to  produce  apparent  temperature  output  files  for  the  plotting  program. 


REFD 

90.000 

2.000 

3. 

1. 

SCEN 

SRCD 

5.000 

0.000 

-1000.000 

0.000 

5.00 

0.000 

5.000 

100.000 

5.0 

METD 

SCN3 

GO 

3.000 

0.125 

330.000 

121. 

27.000 

81. 

60.000 

0. 

1009.728 

4.000 

METD 

GO 

3.000 

300.000 

27.000 

60.000 

1009.728 

4.000 

METD 

GO 

3-000 

270.000 

27.000 

60.000 

1009.728 

4.000 

METD 

GO 

3.000 

240.000 

27.000 

60.000 

1009.728 

4.000 

METD 

DONE 

3.000 

210.000 

27.000 

60.000 

1009.728 

4.000 
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OUTPUT: 


»  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 


FIRE  TYPE:  TRUCK 

X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  M 

( 

X  , 

Y 

z 

) 

OBSERVER 

( 

5.0, 

-1000.0, 

5.0 

) 

TARGET 

( 

5.0, 

100.0, 

5.0 

) 

FIRE 

( 

0.0, 

0.0, 

0.0 

) 

METEOROLOGICAL  DATA  : 

WINDSPEED  3.0  M/S  WIND  DIRECTION  330.0  DEG 

TEMPERATURE  27.0  C  RELATIVE  HUMIDITY  60.0  % 

AIR  DENSITY  1009.7  G/M**3  PASQUILL  CATEGORY  4 


IMAGER  DATA: 

TOTAL  FOV:  15.13  BY  10.13  (MR)  RESOLUTION  PER  PIXEL  0.13  MR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  121,  81  ) 

WAVELENGTH  3.60  UM 


PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.74  M/S  COEFFICIENT  (C)  7.530 

BUOYANCY  FLUX  (F)  41.5  VIRTUAL  SOURCE  POSITION  (H) 

ENTRAINMENT  CONSTANT  0.3820  XV:  -6.17,  ZV:  -8.45 

**************** 

*  FGLOW  OUTPUT  * 

**************** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR  121  COLUMN,  81  ROW  ARRAY. 
THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  -2.50,  0.00,  10.00  )  H 

AND  THE  ARRAY  CENTER  IS  AT  (  5.00,  0.00,  5.00  )  M 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SIXJRCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.125  M  UNITS:  Watts/in**2/sr/um 
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****************** 


*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  M 

( 

X  , 

Y 

z 

) 

OBSERVER 

( 

5.0, 

-1000.0, 

5.0 

) 

TARGET 

( 

5.0, 

100.0, 

5.0 

) 

FIRE 

( 

0.0, 

0.0, 

0.0 

) 

METEOROLOGICAL  DATA  : 

UINDSPEED  3.0  M/S  WIND  DIRECTION  300.0  DEG 

TEMPERATURE  27.0  C  RELATIVE  HUMIDITY  60.0  % 

AIR  DENSITY  1009.7  G/M**3  PASQUILL  CATEGORY  4 

IMAGER  DATA: 

TOTAL  FOV:  15.13  BY  10.13  (MR)  RESOLUTION  PER  PIXEL  0.13  MR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  121,  81  ) 

WAVELENGTH  3.60  UM 


PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.74  M/S  COEFFICIENT  (C)  7.530 

BUOYANCY  FLUX  (F)  41.5  VIRTUAL  SOURCE  POSITION  (M) 

ENTRAINMENT  CONSTANT  0.3820  XV:  -6.17,  ZV:  -8.45 

**************** 

*  FGLOW  OUTPUT  * 

**************** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR  121  COLUMN,  81  ROW  ARRAY. 

THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  -2.50,  0.00,  10.00  )  M 

AND  THE  ARRAY  CENTER  IS  AT  (  5.00,  0.00,  5.00  )  M 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SOURCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.125  M  UNITS:  Watts/ni**2/sr/uni 


****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 

X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  M 

( 

X  , 

Y 

z 

) 

OBSERVER 

( 

5.0, 

-1000.0, 

5.0 

) 

TARGET 

( 

5.0, 

100.0, 

5.0 

) 

FIRE 

< 

0.0, 

0.0, 

0.0 

) 

METEOROLOGICAL  DATA  : 

WINDSPEED  3.0  M/S  WIND  DIRECTION  270.0  DEG 

TEMPERATURE  27.0  C  RELATIVE  HUMIDITY  60.0  % 

AIR  DENSITY  1009.7  G/M**3  PASQUILL  CATEGORY  4 

IMAGER  DATA: 

TOTAL  FOV:  15.13  BY  10.13  (MR)  RESOLUTION  PER  PIXEL  0.13  HR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  121,  81  ) 

WAVELENGTH  3.60  UH 


PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.74  M/S  COEFFICIENT  (C)  7.530 

BUOYANCY  FLUX  (F)  41.5  VIRTUAL  SOURCE  POSITION  (M) 

ENTRAINMENT  CONSTANT  0.3820  XV:  -6.17,  ZV:  -8.45 

**************** 

*  FGLOW  OUTPUT  * 

**************** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR  121  COLUMN,  81  ROW  ARRAY. 

THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  -2.50,  0.00,  10.00  )  M 

AND  THE  ARRAY  CENTER  IS  AT  (  5.00,  0.00,  5.00  )  M 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SOURCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.125  H  UNITS:  Watts/m**2/sr/uni 
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****************** 

*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X-AXIS  HEADING:  90.0  DEG  CW  FROM  NORTH 


POSITION  M 

( 

X  , 

Y 

z 

) 

OBSERVER 

( 

5.0, 

-1000.0, 

5.0 

) 

TARGET 

( 

5.0, 

100.0, 

5.0 

) 

FIRE 

( 

0.0, 

0.0, 

0.0 

) 

METEOROLOGICAL  DATA  : 

WINDSPEED  3.0  M/S  WIND  DIRECTION  240.0  DEG 

TEMPERATURE  27.0  C  RELATIVE  HUMIDITY  60.0  % 

AIR  DENSITY  1009.7  G/M**3  PASQUILL  CATEGORY  4 

IMAGER  DATA: 

TOTAL  FOV:  15.13  BY  10.13  (MR)  RESOLUTION  PER  PIXEL  0.13  MR 
PIXELS  (HORIZONTAL,  VERTICAL)  (  121,  81  ) 

WAVELENGTH  3.60  UM 


PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.74  M/S  COEFFICIENT  (C)  7.530 

BUOYANCY  FLUX  (F)  41.5  VIRTUAL  SOURCE  POSITION  (M) 

ENTRAINMENT  CONSTANT  0.3820  XV:  -6.17,  ZV:  -8.45 

**************** 

*  FGLOW  OUTPUT  * 

**************** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR  121  COLUMN,  81  ROW  ARRAY. 
THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  -2.50,  0.00,  10.00  )  H 

AND  THE  ARRAY  CENTER  IS  AT  (  5.00,  0.00,  5.00  )  M 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SOURCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.125  M  UNITS:  Watts/m**2/sr/uni 


****************** 


*  EOSAEL  92  * 

*  4D  FITTE  INPUT  * 
****************** 


FIRE  TYPE:  TRUCK 


X-AXIS  HEADING: 

90.0  DEG  CW  FROM 

NORTH 

POSITION  M 

( 

X  ,  Y  , 

Z  ) 

OBSERVER 

C 

5.0,  -1000.0, 

5.0  ) 

TARGET 

( 

5.0,  100.0, 

5.0  ) 

FIRE 

( 

0.0,  0.0, 

0.0  ) 

METEOROLOGICAL  DATA  : 

WINDSPEED  3.0  M/S  WIND  DIRECTION  210.0  DEG 

TEMPERATURE  27.0  C  RELATIVE  HUMIDITY  60.0  % 

AIR  DENSITY  1009.7  G/M**3  PASQUILL  CATEGORY  4 

IMAGER  DATA: 

TOTAL  FOV:  15.13  BY  10.13  (MR)  RESOLUTION  PER  PIXEL  0.13  HR 
PIXELS  (HORIZONTAL,  VERTICAL)  (121,  81  ) 

WAVELENGTH  3.60  UM 


PLUME  CONSTANTS  : 

VERTICAL  WIND  (WO)  2.74  H/S  COEFFICIENT  (C)  7.530 

BUOYANCY  FLUX  (F)  41.5  VIRTUAL  SOURCE  POSITION  (M) 

ENTRAINMENT  CONSTANT  0.3820  XV:  -6.17,  ZV:  -8.45 

**************** 

*  FGLOW  OUTPUT  * 

**************** 

AT  100.  SECONDS  AFTER  START  OF  FIRE 

RADIANCE  AND  TRANSMITTANCE  DATA  STORED  FOR  121  COLUMN,  81  ROW  ARRAY. 

THE  TOP  LEFT  CORNER  PIXEL  IS  AT  (  -2.50,  0.00,  10.00  )  M 

AND  THE  ARRAY  CENTER  IS  AT  (  5.00,  0.00,  5.00  )  H 

IN  A  PLANE  (  PERPENDICULAR  TO  THE  LOS  AND  AT  THE  SOURCE  DISTANCE 
FROM  THE  OBSERVER  )  IN  THE  USER  COORDINATE  SYSTEM. 

POINT  SPACING  IS  0.125  M  UNITS:  Watts/in**2/sr/um 


Figure  11.  Apparent  temperature  plots  for  various  wind  directions.  Prepared 
from  the  data  files  of  the  example  9  simulation.  Wavelength  was  3.6  fim. 
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IR 
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Atmospheric  Sciences  Laboratory 

Battlefield  Emissive  Sources  Trial,  Oldebruck,  the  Netherlands, 
Europe 

Battlefield  Induced  Contamination  Test  I 
Battlefield  Induced  Contamination  Test  III 
Electro-Optical  Systems  Atmospheric  Effects  Library 
Fire-Induced  Transmittance  and  Turbulence  Effects 
infi-ared 
line  of  sight 
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Appendix  A 

Fluctuation  Time-History  Data  Files 


A-1.  Basic  Ideas  Underlying  the  Fluctuation  Data  Files 


Fire  plumes  move  about  in  space  as  a  function  of  time.  The  fluctuations  of  the 
centerline  position  at  a  given  plume  height  are  a  result  of  fluctuations  in  the 
horizontal  windspeed  and  wind  direction  since  that  plume  segment  was 
generated.  The  centerline  fluctuations  can  be  described  by  a  deterministic 
model  if  the  wind  fluctuations  are  known.  The  fluctuations  in  temperature  and 
effluent  concentrations  at  any  point  in  the  plume  are  the  result  of  turbulent 
mixing.  For  a  simple  fast-rurming  model,  the  fluctuations  can  be  approximated 
by  a  distribution  with  a  turbulent  power  spectral  density  with  amplitudes  and 
phase  angles  based  on  random  numbers.  The  amplitude  fluctuations  can  be 
scaled  using  available  data. 

The  minimum  information  needed  to  model  a  space-  and  time-varying  (4D) 
plume  is  a  set  of  random  numbers  and  x-  and  y-wind  fluctuation  data. 

Several  additional  considerations  related  to  model  simplicity  were  used  in 
constructing  the  files: 

•  A  1-s  time  step  was  defined  for  the  calculations.  The  row  number  in  the  file 
corresponds  to  the  number  of  seconds  after  the  start  of  the  fire.  The  data  file 
must  be  long  enough  to  describe  plume  transport  after  the  fire  bums  out,  so  a 
total  file  length  of  2000  rows  was  chosen. 

•  It  is  desirable  to  be  able  to  use  the  fluctuation  data  files  for  a  range  of 
windspeeds,  so  the  wind  fluctuation  data  should  be  normalized;  the  data  stored 
should  be  fi'actional  fluctuations. 

•  Calculations  of  centerline  displacement  are  based  on  the  cumulative 
fluctuations  since  the  displaced  plume  segment  was  created.  These  calculations 
are  most  easily  and  quickly  performed  if  the  wind  fluctuation  data  are  stored 
as  cumulative  fi'actional  fluctuations. 

If  model  performance  is  to  be  compared  with  data  from  field  tests,  it  is 
desirable  to  use  a  time-history  file  that  reflects  the  wind  conditions  during  the 
test.  Constmction  of  such  a  file  is  described  below. 
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A-2.  Construction  of  a  Time-History  File  from  Field  Test  Data 


A  time-history  file  may  be  constructed  from  field  test  data  if  data  meet  the 
following  criteria: 

•  include  windspeed  and  wind  direction 

•  available  at  1-s  intervals 

•  contain  (at  least)  2000  observations. 

The  2000  s  of  data  must  be  averaged  to  obtain  mean  windspeed  and  wind 
direction.  The  mean  windspeed  and  wind  direction  should  be  recorded  for  use 
as  input  on  the  METD  input  card.  The  data  must  be  processed  to  obtain  time- 
histories  of  windspeed  in  and  perpendicular  to  the  mean  wind  direction. 

The  data  can  be  transformed  and  stored  in  a  time-history  file.  The  file,  as 
mentioned  earlier  in  the  manual,  must  have  three  columns  of  data.  The  first 
column  must  contain  a  set  of  random  numbers  in  the  range  -0.5  to  0.5,  and  a 
mean  of  0.00  +  0.009.  The  second  column  must  contain  the  cumulative 
normalized  fluctuations  in  the  downwind  component  of  the  wind.  If  we  define 
the  mean  windspeed  as  U  and  the  instantaneous  windspeed  in  the  mean  wind 
direction  as  the  value  in  the  n*  row  of  the  second  column  is  given  by 


The  third  column  of  the  data  file  contains  the  cumulative  normalized  crosswind 
component.  If  we  define  the  instantaneous  value  of  the  crosswind  component 
to  be  v„  the  value  in  the  n*  row  of  the  third  column  is  given  by 


(A-2) 


The  time-history  file  is  read  with  an  unformatted  READ  statement  that  requires 
three  values  for  each  row,  so  the  columnar  values  for  each  row  should  be 
separated  by  spaces  or  other  separators  recognized  by  FORTRAN. 
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You  should  rename  the  existing  file  FIT4D.DAT  and  copy  the  desired  time- 
history  file  into  FIT4D.DAT,  because  it  is  the  filename  used  by  FITTE  for  the 
time-history  file. 


Appendix  B 

Auxiliary  Programs  for  Visualization 
of  FGLOW  Output 


Two  accessory  programs  are  provided  to  help  with  visualization  of  the  image 
files  created  by  FITTE.  The  accessory  programs  are  working  examples  that 
transform  the  FGLOWn.FDT  files  into  forms  that  may  be  printed  or  imported 
into  a  graphics  program  that  can  produce  contour  plots.  Both  of  the  programs 
permit  output  of  four  quantities:  radiance,  transmittance,  apparent  (blackbody- 
equivalent)  temperature,  and  integrated  path  density  (CL).  The  latter  two 
quantities  are  obtained  by  transformation  of  the  first  two  for  the  appropriate 
wavelength.  (This  would  be  the  central  wavelength  for  a  waveband 
calculation.) 

SEPARATE  SIMULATION  RUNS  ARE  NO  LONGER  NECESSARY  TO 
OBTAIN  APPARENT  TEMPERATURE  IMAGES. 

The  FGLOWn.FDT  output  files,  generally,  should  be  renamed  for  storage; 
otherwise,  they  will  be  overwritten  the  next  time  FITTE  is  run  using  the 
FGLOW  identifier. 

The  programs  are  written  for  use  with  equipment  and  software  that  we  have  on 
hand.  This  should  not  be  construed  as  a  recommendation  or  endorsement  of  the 
hardware  and  software. 

Plots  prepared  with  output  from  XS2GRD  were  shown  in  section  5  with  sample 
runs.  Appendix  B  contains  a  sample  of  output  from  XS2PRT. 

B-L  XS2PRT 

XS2PRT  creates  files  containing  greyscale  printouts  of  the  FGLOWn.FDT  data. 
The  files  have  the  same  name  as  the  input  file  with  the  extension  FDT  replaced 
by  RAD,  TRN,  APT,  or  DNS  as  appropriate.  The  example  program  is  tailored 
for  an  HP  LaserJet  printer  and  uses  control  sequences  and  built-in  fonts  in 
creating  the  output  files.  The  width  of  the  image  array  (that  is,  the  number  of 
columns  in  the  array)  determines  whether  the  output  is  printed  in  portrait  or 
landscape  mode  or  is  split  into  panels  for  printing. 

If  another  printer  is  used  for  the  printouts,  it  is  necessary  to  modify  the 
program.  This  entails  substitution  of  the  appropriate  printer  control  sequences 
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(usually  found  in  the  printer  manual)  and  modification  of  the  numbers  that 
determine  when  to  split  the  output  into  panels.  Fixed  spacing  fonts  should  be 
used,  and  the  number  of  characters  per  inch  should  be  (nearly)  equal  to  the 
number  of  lines  per  inch  to  give  die  image  the  correct  proportions.  The 
program  is  relatively  well  commented  to  assist  in  identifying  areas  needing 
changes. 

Figure  B-1  shows  a  printout  of  integrated  path  density  created  using  this 
program.  The  simulation  file  was  the  same  file  used  to  create  the  plume 
transmittance  contour  image  shown  in  figure  4.  The  greyscale  assignment  table 
was  moved  to  the  same  page  as  the  image  for  this  example. 

B-2.  XS2GRD 

XS2GRD  transforms  the  FGLOWn.FDT  files  to  a  form  that  commercial 
software  uses  to  create  surface  or  contour  plots.  This  program  should  be 
relatively  easy  to  modify  for  use  with  other  plotting  packages.  The  format  for 
the  output  file  is  described  in  program  comments. 

Because  the  program  used  for  the  contour  plots  expects  a  file  extension  of 
.GRD,  the  output  files  for  this  program  have  the  extension  .FDT  replaced  by 
R.GRD,  T.GRD,  A. GRD,  or  D.GRD  depending  on  which  of  the  four  types  of 
data  it  contains. 
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SIKULATIOK  FILE  IS  C:\FIG3M.FDT  :  PLUME  DEKSITY  (CL) 
0.630000  uM,  IMAGE  WIDTH,  HEIGHT  (M) :  60.0000 


FOR  WAVELENGTH 
,  40.0000 


KEY  TO  IMAGE  VALUES: 
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0.390212  TO  0.445957 
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0.613190  TO  0.668935 

0.668935  TO  0.724680 

0.724680  TO  0.780424 

0.780424  TO  0.836169 

0.836169  TO  0.891913 


Figure  B-1.  Greyscale  print  of  plume  optical  density  produced  with  XS2PRT. 
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Appendix  C 

Modification  of  Fire  Parameters 


The  selection  of  fire  parameters  to  provide  realistic  representations  within  the 
model  is  not  a  simple,  straightforward  process.  There  are  also  limitations, 
primarily  imposed  by  model  assumptions  that  limit  the  fires  that  can  be 
described.  Appendix  C  discusses  the  limitations,  describes  rough  relations 
between  model  and  physical  fire  parameters,  and  suggests  limits  on  parameter 
variation  using  the  SVAR  input  card. 

C-1.  Model  Limitations  on  Fire  Parameters 

FITTE  was  developed  as  a  fast  running  code  to  model  effects  of  small-scale 
battlefield  fires  that  were  assumed  to  be  vehicular  fires.  Table  C-1  shows  basic 
assumptions  made  (consistent  with  this  objective)  and  associated  model 
limitations.  The  first  and  second  assumptions  must  be  considered  when  fire 
parameters  are  varied.  The  size  limits  shown  are  approximate.  It  is  difficult 
to  set  bounds  for  the  heat  output,  but  the  values  associated  with  the  fires 
defined  in  the  model  are  the  right  order  of  magnitude. 


Table  C-1.  Limitations  caused  by  model  assumptions 


Model  Assumption 

Associated  Limitations 

1 .  The  energy  density  of  the  fire  plume 
is  low  compared  with  that  of  the 
atmosphere. 

Only  small-scale  fires  (  R  <  10  m) 
with  low  to  moderate  heat  output  can 
be  handled. 

2.  The  plume  must  have  structural 
integrity. 

The  fire  has  a  minimum  size 
(  R  >  0.5  m  )  and  heat  flux. 

3.  The  mean  plume  centerline  is  described 
by  a  power  law. 

Unusual  wind-shear  conditions 
cannot  be  handled. 

Less  accuracy  under  very  low 
windspeed  conditions. 

4.  Fires  are  assumed  to  be  circular. 

Cannot  easily  be  adapted  to  describe 
linear  fires,  such  as  grass  fires. 
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The  limits  are  intended  to  indicate  that  model  predictions  have  built-in  unknown 
errors  if  parameters  are  varied  to  describe  other  fires,  like  an  oil  well  fire 
because  it  is  a  high-energy  density  source. 

The  parameters  that  can  be  varied  using  the  input  card  SVAR  are  mean  fire 
temperature,  aerosol  efficiency  factor,  fire  radius,  and  bum  time.  It  is  possible 
to  vary  these  parameters  to  examine  the  effect  of  uncertainties  associated  with 
their  values  or  to  provide  a  better  comparison  with  field  test  data. 

Variation  of  the  parameters  to  examine  the  effect  of  uncertainties  associated 
with  their  values  was  discussed  briefly  in  section  4.  For  providing  a  better 
comparison  with  field  test  data,  the  following  comments  may  be  useful. 


C-1.1  Developing  Model  Fire  Parameters  from  Data 

The  bum  time  can  be  obtained  fi’om  observations.  The  associated  uncertainty 
lies  in  determining  the  end  time  for  the  fire.  The  end  time  should  be  chosen 
as  the  time  when  a  cohesive  plume  ceases  to  exist. 

The  fire  radius  has  a  clear  meaning  for  a  pool  fire  in  a  circular  tank.  If  another 
container  shape  is  involved,  such  as  a  rectangular  tank  of  dimensions  L  by  W 
with  L  <  2W,  an  appropriate  value  would  probably  be  R  =  0.5  (L  +  W).  The 
exact  shape  is  not  critical  because  the  influx  of  air  modifies  the  shape  of  the 
fire  or  fire  plume  so  it  is  quasicircular. 

The  mean  fire  temperature  can  be  estimated  firom  a  set  of  measurements  above 
the  fire.  This  can  be  tricky  for  a  fire  in  the  field,  and  the  estimate  will  have  a 
relatively  large  uncertainty.  The  value  used  in  the  model  was  obtained  as 
follows  (using  Kelvin  temperatures): 

1.  The  data  were  examined  to  obtain  the  (average  of  the)  maximum  observed 
temperatures  for  the  thermocouples  near  the  center  of  the  fire.  (The 
thermocouple  array  was  at  about  the  level  of  the  flame  tips.)  The  (average) 
value  was  assumed  to  be  valid  for  a  point  at  a  radius  of  la  fi-om  the  center  of 


101 


a  Gaussian  distribution.  (The  wind  shifted  the  fire  about,  and  the  thermocouple 
array  was  in  an  area  with  a  narrow  distribution  function.) 

2.  The  value  was  divided  by  exponent(-.5)  to  obtain  an  estimate  of  the 
Gaussian  maximum. 

3.  The  ambient  temperature  was  subtracted  from  the  Gaussian  maximum 
temperature  to  obtain  the  excess  plume  maximum  temperature.  The  excess 
plume  maximum  temperature  was  divided  by  2.9;  the  ratio  of  Gaussian/top  hat 
distribution  values  for  the  observation  conditions.  The  ambient  temperature  was 
then  added  to  obtain  the  mean  fire  temperature. 

It  is  extremely  difficult  to  obtain  a  good  measurement  of  the  aerosol  efficiency 
factor.  The  value  built  into  FITTE  should  be  appropriate  under  most 
conditions. 

C-1.2  Checking  Model  Behavior  When  Varying  Fire  Parameters 

If  plume  image  data  (photographs,  video,  or  infrared  images)  are  available,  it 
is  possible  to  check  the  plume  behavior  to  ensure  that  the  fire  parameters 
produce  a  realistic  plume.  The  following  checks  may  be  made  by  comparing 
the  data  with  information  derived  from  the  model  output.  The  most  accurate 
checks  are  based  on  images  taken  from  a  crosswind  position. 

The  quantities  to  be  checked  are  the  initial  plume  slope  and  spread.  A 
scenario  1  model  run  with  the  fire  parameter  variation  produces  output  that  may 
be  used  to  predict  these  values. 

The  ratio  of  the  vertical  wind  to  the  windspeed  gives  the  initial  plume  slope. 
The  inverse  tangent  predicts  the  initial  elevation  angle.  The  elevation  angle 
should  be  compared  with  the  slope  or  elevation  angle  obtained  from  the  plume 
images.  A  short  segment  of  the  plume  image  should  give  the  best  value  for  the 
angle.  (You  may  want  to  try  this  comparison  with  example  10  to  get  an  idea 
of  the  errors  involved.) 
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The  plume  spread  angle  is  the  inverse  tangent  of  the  entrainment  constant.  You 
should  compare  the  plume  spread  angle  with  the  (full  width)  angular  spread  of 
the  plume  obtained  from  the  images.  Use  as  long  a  plume  segment  as  possible 
for  the  most  accurate  result. 

If  agreement  is  good  for  both  comparisons,  the  model  should  provide  good 
predictions  of  the  propagation  effects  when  the  selected  fire  parameters  are 
used.  Good  agreement  for  the  initial  elevation  angle  depends  on  the  variability 
of  the  wind  during  the  field  test  as  well  as  the  uncertainty  in  determining  the 
angle  from  the  images.  The  plume  spread  angle  should  agree  within  5°. 

A  scenario  3  run  could  also  be  used;  a  printed  simulated  image  could  be 
produced  with  the  accessory  program  XS2PRT  and  could  be  compared  directly 
with  the  real  plume  images. 
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Appendix  D 
Model  Organization 


Figure  D-1  shows  the  organization  of  FITTE.  It  can  be  called  from  EOSAEL 
main,  or  run  as  a  stand-alone  driver  called  by  FITDRV .  FOR.  If  FITDRV .  FOR 
is  used,  it  and  its  corresponding  block  data  subroutine  BLOOD .  FOR  must  be 
customized  for  the  system  on  which  it  is  to  be  run. 

In  general,  the  subroutines  called  directly  by  FITTE  are  called  once,  except  for 
XSMIT  which  is  called  once  for  each  path  segment.  For  scenario  3,  SETLOS 
is  called  to  establish  parameter  values  for  each  line  of  sight  (LOS)  of  the  image 
matrix. 

The  FITTE  model  consists  of  26  subroutines  and  12  functions.  EOSAEL  main 
or  FITDRV.  FOR  is  used  as  the  main  program,  and  the  corresponding  block 
data  routines,  FITTED  and  BLOOD,  respectively,  must  be  linked  with  the 
program. 

A  brief  description  of  each  subroutine  and  function  is  given  below.  The 
subroutines  are  listed  roughly  as  they  are  first  called.  Subroutines  that  perform 
specific  related  calculations  are  grouped  together. 

D-1.  Subroutines 

D-l.l  FITTE 

FITTE  is  the  subroutine  that  controls  execution  of  the  model.  It  calls 
subroutines  to  perform  required  calculations  and  write  model  output  files. 

D-1. 2  FITTED 

FITTRD  controls  reading  of  the  FITTE  input  cards,  checks  the  input  values, 
substitutes  defaults  if  possible  when  invalid  data  are  encountered,  and  calls 
subroutines  to  transform  coordinates  and  units  and  calculate  plume  constants  for 
given  input  parameters. 

D-1. 2. 1  XSFORM. — ^XSFORM  transforms  EOSAEL  common  block  input  to  units  used  by 

FITTE. 
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Figure  D-1.  Diagram  of  FITTE  model  structure.  Dashed  lines  indicate  conditional 
calls.  Functions  called,  except  EXPOF,  are  listed  in  parentheses  in  subroutine  boxes. 

D-1. 2. 2  ROTOR. — ROTOR  transforms  coordinates  to  (and  from)  the  wind-aligned  system 

used  by  FITTE. 

D-1. 2. 3  SOURC. — SOURC  sets  values  of  fire  parameters  based  on  input  from  the  REFD 

card  and  the  SVAR  card.  Table  1  of  this  document  shows  the  parameters  for 
the  various  sources. 

D-1. 2. 4  SBAR2. — SBAR2  calculates  the  stability  parameter  needed  to  calculate  the 

ambient  lapse  rate  for  FITTE.  The  routine  is  identical  with  subroutine  SBARl 
used  in  the  EOSAEL  82  version  of  model  COMBIC. 


D-1. 2. 5  FITPLM. — FITPLM  calculates  the  values  of  FITTE  plume  parameters. 
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D-1.3 


GLOWA 


D-1.4 

D- 1.4.1 

D-1.4. 2 

D-1.5 

D-1.6 

D-1.7 

D-1.7.1 


GLOWA  is  called  to  calculate  parameters  needed  for  control  of  FGLOW 
(ISCN  =  3)  model  execution. 

CHITMP 

CHITMP  calculates  values  of  plume  distribution  functions,  such  as  temperature 
and  aerosol  and  molecular  densities  for  a  given  point  on  the  LOS. 

FLUCTl. — FLUCTl  calculates  the  displacement,  caused  by  wind  fluctuations, 
of  a  plume  centerline  point  from  the  mean  centerline  for  the  current  calculation 
time. 

FLUCT2. — FLUCT2  calculates  fluctuations  in  entrainment  and  combustion 
products  caused  by  wind  fluctuations. 

SETLOS 

SETLOS  calculates  plume  entry  and  exit  points  on  the  LOS. 

XSMIT 

XSMIT  controls  calculation  of  the  propagation  parameters  related  to  extinction 
and  radiance.  When  molecular  absorption  effects  are  to  be  included 
(IMOL  o  0),  it  calls  the  set  of  subroutines  adapted  from  the  ATLES  model 
code. 

ATLESF 

ATLES F  controls  calculations  of  molecular  absorption  by  the  combustion 
products  carbon  dioxide  and  water  vapor.  Temperature-dependent  absorption 
coefficients  are  used  for  the  calculations.  Values  are  calculated  for  each  of  the 
NWYN  wavenumbers  if  a  waveband  calculation  is  specified. 

BAND. — BAND  reads  the  temperature-dependent  band  model  parameters. 
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D-1.7.2 

D-1.7.3 

D-1.7.4 

D-1.8 

D-L9 

D-1.10 

D-Lll 

D- 1.1 1.1 


SPCLOP. — SPCLOP  creates  an  array  IWAVE  to  control  the  band  model 
calculations. 

BP  ARAM. — ^BPARAM  calculates  band  model  parameters  for  a  given  wave 
number,  path  segment,  and  species  of  molecule. 

DLTRAN. — ^DLTRAN  calculates  the  transmittance  along  a  path  segment  for  a 
single  molecular  species. 

TRAN 

TRAN  averages  transmittances  and  integrates  radiances  over  the  waveband.  It 
calls  RESPNS  to  create  an  instrument  response  function  array  for  use  in  these 
calculations,  and  it  calls  PRCOEF  to  obtain  the  aerosol  optical  constants  for 
each  wavelength.  For  single  wavelength  calculations,  it  transfers  values  from 
the  calculation  arrays  to  the  output  variables. 

RESPNS 

RESPNS  calculates  an  array  of  instrument  response  coefficients  for  the  wave 
numbers  used  in  a  waveband  calculation. 

PRCOEF 

PRCOEF  assigns  values  of  aerosol  optical  constants  for  a  given  wavelength. 
FIJIT 

FIJIT  controls  calculation  of  the  optical  turbulence  effects  parameters.  It  calls 
one  or  both  of  the  routines  listed  below,  depending  on  the  FITTE  scenario 
specified. 

COHERE. — COHERE  calculates  the  lateral  coherence  length  by  using  ftmction 
ASRI  to  integrate  the  function  FRHOO  through  the  plume. 
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D- 1.1 1.2 

D-L12 

D- 1.1 2.1 

D-2. 

D-2.1 

D-2.2 

D-2.3 


RADII. — RADII  performs  the  calculations  described  by  Fante  (1975),  using 
ASRI  to  perform  the  necessary  integrations.  It  assigns  a  value  to  I  CASE  to 
indicate  which  of  the  four  calculation  methods  was  used. 

FITTWT 

FITTWT  controls  formatting  and  printing  of  user  input  and  model  output.  The 
output  is  tailored  to  the  specified  scenario.  A  warning  is  printed  if  defaults 
were  substituted  for  any  of  the  user-specified  variables,  and  an  error  message 
is  printed  if  fatal  input  data  are  detected.  The  value  of  MODE  is  used  to  specify 
whether  input  or  output  is  written. 

PRTSPC. — PRTSPC  prints  a  table  of  spectrally  resolved  transmittance  and 
radiance  data. 

Functions 

APPTMP 

APPTMP  is  an  inverse  blackbody  function  used  to  produce  apparent  temperature 
output  for  FGLOW  if  I  APT  =1.  For  a  waveband  calculation,  the  evaluation 
is  performed  for  the  wavelength  at  the  center  of  the  band. 

ASRI 

ASRI  is  an  adaptive  Simpson’s  rule  integration  routine  with  a  specified  local 
error  and  a  maximum  permitted  number  of  resolution  increases.  If  the  function 
being  integrated  does  not  meet  the  error  tolerance  after  the  maximum  number 
of  resolutions  increases  for  a  given  path  interval,  a  flag  is  set  to  indicate  that 
the  integration  error  is  greater  than  that  specified. 

BLACKS 

BLACKB  calculates  the  energy  density  of  a  blackbody  of  temperature  T  at 
wavelength  X. 
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D-2.4 


EXPOF 


EXPOF  limits  the  allowed  range  of  exponents  to  prevent  program  crashes 
caused  by  calculation  errors. 

D-2.5  FCN2 

FCN2  evaluates  C/  as  a  function  of  temperature  at  any  point  in  the  plume. 
D-2.6  FLUCT3 

FLUCT3  evaluates  the  local  fractional  fluctuation  in  plume  temperature  based 
on  the  assumption  of  a  fractal  spatial  frequency  distribution. 

D-2.7  FLUCT4 

FLUCT4  evaluates  the  local  fractional  fluctuation  in  plume  effluents  based  on 
the  assumption  of  a  fractal  spatial  frequency  distribution. 

D-2.8  FRHOO 

FRHOO  evaluates  the  integrand  in  the  equation  for  RHOO  in  subroutine  COHER. 
D-2.9  FRHOL 

FRHOL  evaluates  the  integrand  in  the  equation  for  Pl  in  subroutine  RADII. 

D-2.10  MU 

MU  evaluates  the  ratio  of  short-  to  long-term  averaged  beam  radius  using 
numerical  results  (Fante  1975).  This  function  is  used  by  subroutine  RADII  for 
case  2  calculations  as  described  in  the  technical  manual. 
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D-2.11  QARTIC 

QARTIC  solves  a  quartic  equation  to  find  the  eenterline  point  closest  to  a  given 
LOS  point.  It  is  called  by  CHITMP. 

D-2J2  RTVAR 

RTVAR  evaluates  the  variance  of  1/T  for  function  FCN2  of  subroutine  RADII. 
It  uses  a  look-up  table  and  an  interpolation  scheme. 
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